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ABSTRACT 

A  high-resolution,  multi-level,  primitive  equation  ocean  model  is  used  to  examine 
the  response  of  an  idealized,  flat-bottom,  oceanic  regime  off  northern  California  to 
steady,  equatorward,  local  wind-forcing  during  the  upwelling  season.  The  model  has 
open  boundaries  on  all  but  the  eastern  coastal  boundary  on  which  either  free-slip  or 
zero-slip  boundary  conditions  are  imposed.  Time-invariant  winds,  either  with  or 
without  a  component  of  wind  stress  curl,  are  used  as  model  forcing  to  spin-up  a 
classical  two-dimensional,  upwelling-induced  coastal  jet  and  undercurrent.  Since  no 
eddies  are  generated,  a  stability  analysis  of  the  mean  flow  is  conducted  which  explores 
both  the  necessary  conditions  of  mixed  (barotropic  and  baroclinic)  instability,  through 
calculations  of  potential  vorticity,  and  the  sufficient  conditions  for  baroclinic 
instability,  through  an  application  of  a  simple  two-layer  stability  model.  Comparisons 
of  model  results  with  observations  of  the  coastal  jet  in  the  California  Current  System 
indicate  that  the  location  and  the  horizontal  and  vertical  current  shear  associated  with 
the  model  coastal  jet  compare  favorably  with  observations;  however  the  modeled  jet  is 
stronger,  deeper  and  wider  than  the  observed  jet.  Finally,  the  inclusion  of  wind  stress 
curl  and  the  zero-slip  boundary  condition  are  demonstrated  to  be  important  elements 
in  model  simulations  of  the  coastal  jet. 
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I.  INTRODUCTION 

A.       OBJECTIVES 

The  attainment  of  sufficient  knowledge  to  adequately  simulate  observed  synoptic- 
mesoscale  features  of  the  California  Current  System  (CCS),  and  to  understand  the 
intricate  processes  responsible  for  spatial  and  temporal  variations  within  the  system,  is 
the  goal  of  a  numerical  modeling  effort  presently  underway  at  the  Naval  Postgraduate 
School  in  Monterey,  California.  The  ability  to  predict  positions,  intensities,  circulation 
patterns,  etc.  of  fronts,  undercurrents,  eddies  and  current  filaments,  given  an  initial  set 
of  observations,  could  provide  invaluable  assistance  to  modern  naval  operations 
concerned  with  boundary  current  acoustic  characteristics  and  associated  naval  tactics. 

Recently,  U.S.  Naval  battle  group  tactics  have  been  modified  by  the  predictions 
of  a  quasigeostrophic  (QG)  model  developed  at  Harvard  University  in  conjunction  with 
the  Naval  Postgraduate  School.  The  QG  model  is  initialized  by  infra-red  satellite 
imagery  and  air  dropped  bathythermograph  transects  across  the  Gulf  Stream  and 
prominent  Gulf  Stream  meanders  and  rings.  After  initialization,  the  QG  model  is  able 
to  predict  oceanic  conditions  out  to  approximately  one  week.  Accuracy  is  further 
enhanced  by  daily  model  updates  and  m-situ  verification  procedures.  The  initial  results 
of  this  operation  look  extremely  promising  and  there  is  potential  for  expansion  to  other 
operating  areas  of  the  world.  With  the  ever  increasing  capabilities  of  high  speed 
computers,  the  possibility  of  using  a  high  resolution  Primitive  Equation  (PE)  model  in 
an  eastern  boundary  current  regime,  for  instance,  is  feasible.  The  vorticity  signature  of 
an  eastern  boundary  current  is  generally  less  than  that  of  a  western  boundary  current 
due  to  the  absence  of  western  intensification,  and  the  observed  scale  of  eddies  is  also 
smaller.  Additionally,  while  the  Gulf  Stream,  especially  its  extension  region,  is  in  part 
an  open  ocean  jet  and  frontal  system,  the  CCS  is  very  much  associated  with  the  coastal 
upwelling  region.  This  region  is  characterized  by  significant  horizontal  variations  in 
stratification  and  large  (topographically  and  coastally  induced)  vertical  velocities.  For 
these  reasons,  a  high  resolution  PE  model  may  be  more  applicable  than  a  QG  model 
for  use  in  areas  such  as  the  California  Current  region.  Other  agencies,  as  well,  could 
greatly  benefit  from  such  a  predictive  model,  such  as  the  Coast  Guard  in  search  and 
rescue  efforts  or  agencies  involved  in  oil  spill  recovery,  ocean  engineering,  resource 
exploration,  and  at-sea  waste  disposal  operations. 
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The  intent  of  this  research  is  to  provide  a  contribution  towards  understanding  the 
processes  responsible  for  synoptic-mesoscale  variability  in  the  CCS.  The  tool  used  in 
this  research  effort  is  a  ten-level,  high  resolution  PE  ocean  model  with  idealized 
topographic  capability.  The  study  region  extends  approximately  500  km  offshore  from 
the  west  coast  of  North  America,  and  it  spans  the  California  coastline  from  Point  Sur 
in  the  south  to  Cape  Blanco  in  the  north  (Figure  1.1).  The  particular  realm  of  interest 
is  the  effect  of  wind  forcing  in  the  CCS  regime.  Previous  to  this  effort,  an  imposed 
baroclimc  jet.  intending  to  simulate  the  mean  CCS.  was  used  as  the  forcing  mechanism 
to  explore  potential  instabilities  of  the  mean  flow  (Batteen  et  al,  1985).  Instability  of 
the  mean  flow  occurred  in  both  Hat  bottom  and  topographic  experiments.  It  is 
believed  that  wind  forcing  of  the  local  ocean  domain  could  be  an  equally  important 
process  in  the  CCS  and  may  provide  further  insight  into  the  dynamics,  kinematics,  and 
energetics  of  this  complex  eastern  boundary  current. 

B.       MOTIVATION  AND  BACKGROUND 

The  question  as  to  how  to  design  a  limited  series  of  experiments  to  study  the 
complex  mechanisms  governing  the  eastern  boundary  current  regime  is  not  a  trivial 
one.  Implementation  of  wind  forcing  in  a  numerical  model  must  be  accomplished  m  a 
logical  manner.  The  exact  method  and  design  of  wind  forcing  experiments  raises  many 
questions.  Scientific  foresight  and  results  from  past  wind  forcing  studies  must  be 
intelligently  used,  especially  when  confronted  with  limited  computer  resources. 
Procedurally,  the  most  important  criteria  in  the  design  of  wind  forcing  experiments  is 
probably  to  include  only  one  potential  generating  mechanism  at  a  time,  making  it 
possible  to  isolate  that  effect  during  the  analysis  and  interpretation  phase. 

The  elimination  of  bottom  topography  in  initial  wind  forcing  experiments  allows 
the  role  of  wind  forcing  to  be  isolated  from  the  possible  coupled  role  of  wind  forcing 
with  bottom  topography.  The  particular  developmental  stage  of  the  PE  model  used  in 
this  study  is  such  that  it  does  not  presently  include  the  depth-averaged,  or  barotropic 
flow,  component  (although  it  is  in  the  process  of  being  implemented).  It  is  believed 
that  topographic  effects  would  not  be  well  represented  unless  the  depth-averaged 
component  is  included.  This  does  not.  however,  dismiss  the  effect  of  bottom 
topography  as  an  important  mechanism  in  the  CCS.  The  study  of  bottom  topography 
should  be  given  great  consideration  in  future  development  of  this  numerical  model. 
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It  is  known  that  local  and  remote  winds  drive  coastal  circulation,  the  latter  by 
exciting  coastal  waves  {Carton.  1984).  Wind-driven  coastal  upwelling  ofT  California  is 
closely  tied  to  the  seasonal  cycle  because  northerly  winds  along  the  coast,  which  are 
favorable  for  upwelling,  are  strongest  in  a  climatological  sense  during  summer  (Hickey. 
1979).  Due  to  the  strength  of  the  local  winds  and  the  correlation  of  local  winds  with 
upwelling,  this  study  will  concentrate  solely  on  the  effects  of  local  wind  forcing. 

Chelton  (1984)  points  out  that  the  detailed  temporal  and  spatial  characteristics 
of  the  wind  field  over  the  California  Current  are  not  well  known.  Probably  the  most 
complete  description  of  seasonal  winds  in  the  California  coastal  region  is  given  by 
Nelson  (1977).  who  compiled  direct  ship  observations  in  1°  square  areas  from  records 
dating  back  to  the  mid-nineteenth  century.  Chelton  (19S4)  additionally  notes  that  the 
spatial  structure  of  the  wind  field  is  similar  throughout  the  year  but  that  the  magnitude 
changes  seasonally.  Although  the  time  dependent  nature  of  winds  is  an  important 
consideration,  it  is  seen  here  that  an  important  scale  of  wind  forcing  in  the  CCS  region 
is  on  the  order  of  a  season.  Numerical  studies  to  date  have  considered  many  different 
time  scales  of  forcing.  Carton  and  Philander  (19S4)  state  that  "realistic  coastal  winds 
are  an  essentially  random  succession  of  storms  with  differing  scales."  Philander  and 
Yoon  (1982)  consider  periodic  winds  and  Carton  (1984)  describes  the  coastal  response 
to  an  isolated  storm.  This  study  will  focus  on  the  response  to  forcing  at  the  longest 
time  scale  considered  to  be  of  importance;  that  being  the  seasonal  time  scale.  Shorter 
time  scales  are,  indeed,  interesting  and  warrant  investigation  as  well.  This  study  will 
only  cover  the  summer  season  due  to  the  interest  in  upwelling  mechanisms. 
Furthermore,  because  the  generation  period  of  eddies  (less  than  a  month)  observed  in 
the  CCS  is  short  compared  to  the  seasonal  time  scale,  the  winds  will  be  input  as  time 
invariant  (steady). 

Carton  and  Philander  (1984)  discuss  that  at  seasonal  time  scales,  two 
characteristics  of  the  wind  field  have  been  considered  as  possibly  driving  the  seasonal 
circulation  of  the  eastern  ocean:  the  alongshore  component  of  wind  stress,  and  wind 
stress  curl.  Bryan  and  Ripa  (1978),  Hickey  (1979),  Chelton  (1984).  and  others  have 
suggested  that  upwelling  due  to  the  curl  o[  the  wind  stress  is  an  important  mechanism 
in  establishing  pressure  gradients.  Carton  and  Philander  (1984)  conclude  that  the  wind 
stress  curl  is  important  as  a  forcing  function  and  that  it  significantly  contributes  to 
variability  near  the  coast.  Anderson  and  Gill  (1975)  and  Philander  and  Yoon  (19S2) 
conclude  that  the  alongshore  component  of  the  wind  stress,  which  need  not  contain 
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any  curl,  can  also  contribute  to  the  seasonal  variability  of  the  CCS.  Chelton  (1984) 
points  out  that  the  restriction  to  a  zonally  uniform  wind  stress  could  be  an  important 
weakness  of  eastern  boundary  current  models  because  there  are  strong  cross-shore 
gradients  in  the  wind  stress  resulting  in  a  nearshore  positive  wind  stress  curl.  It  is  seen 
in  the  climatological  wind  record  that  the  winds  over  the  region  of  the  California 
Current  appear  to  be  strongest  approximately  200  km  offshore  so  that  there  is  a 
nearshore  positive  wind  stress  curl  year  round.  Additionally,  it  is  noted  that  the 
magnitude  of  this  wind  stress  curl  vanes  seasonally,  approximately  in  phase  with  the 
alongshore  wind  stress.  In  a  recent  numerical  study  on  the  dynamics  of  the  CCS. 
McCreary  et  ai,  (19S6)  sought  to  answer  the  question  of  what  mechanisms  force 
equatorward  flow  off  the  California  coast  in  a  region  of  positive  wind  stress  curl.  An 
important  result  of  that  study  is  that  not  only  the  alongshore  component  of  the  wind  is 
responsible  for  the  coastal  jet  but  also  that  the  positive  wind  stress  curl  is  another 
source  of  equatorward  flow.  Obviously,  wind  stress  curl  plays  a  crucial  role  in  the 
dynamics  of  the  CCS  and  this  study  will  investigate  that  effect  to  some  extent. 

Another  important  question  which  needs  to  be  addressed  is  whether  the  solutions 
are  sensitive  to  the  frictional  boundary  condition,  i.e.,  free-slip  or  zero-slip,  along  the 
rigid  eastern  boundary  of  the  numerical  model.  From  classical  fluid  dynamics  theory  it 
is  expected  that  the  water  particles  in  direct  contact  with  a  rigid  wall  will  have  zero 
tangential  velocity.  This  zero-slip  condition  has  greater  implications  as  the  resolution 
of  a  particular  model  increases.  The  PE  model  used  in  this  study  has  a  zonal 
resolution  of  8  km  and  the  zero-slip  condition,  therefore,  is  potentially  an  important 
consideration.  This  is  because  a  uniform  long-shore  wind  stress  having  no  curl  will 
produce  vorticity  in  the  ocean  near  the  wall  if  the  condition  of  zero-slip  is  used,  but  no 
vorticity  will  be  produced  if  the  free-slip  condition  is  used. 

In  a  study  of  boundary  conditions  in  homogeneous  ocean  models,  Blandford  (1971) 
explored  the  effect  of  various  boundary  conditions  (free-slip  and  zero-slip)  in  strongly 
nonlinear  flow  in  a  barotropic,  P-plane  model  of  the  ocean.  He  found  that  the  solution 
changes  from  a  steady  flow  for  free-slip  to  an  unsteady  flow  with  the  generation  of 
eddies  along  the  western  boundary-  when  zero-slip  boundary  conditions  are  imposed. 
This  study  will  compare  the  effects  of  free-slip  versus  zero-slip  boundary  conditions  on 
the  eastern  boundarv. 
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With  the  guidance  of  important  results  of  past  studies,  this  study  of  wind  forcing 
of  the  CCS  will  focus  on  local,  steady  wind  forcing  of  a  flat  bottom  regime  with  the 
additional  effect  of  wind  stress  curl.  Also,  both  free-slip  and  zero-slip  boundary 
conditions  will  be  explored.  The  experiments  have  been  designed  to  proceed  in  a 
scientifically  logical  manner  with  the  purpose  of  investigating  one  physical  mechanism 
at  a  time. 
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II.  MODEL  DESCRIPTION 

A.       BASIC  MODEL  DESCRIPTION 
1.  Model  Equations 

The  numerical  model  used  in  this  study  is  a  ten-level  PE  model  of  a  baroclinic 
ocean  on  an  f-plane  with  no  depth-averaged  flow  (the  depth-averaged  component  is 
presently  being  incorporated).  The  model  is  based  on  the  hydrostatic  and  Boussinesq 
approximations  with  the  rigid  lid  approximation  made.  Depth  of  the  model  is  variable, 
when  idealized  bottom  topography  is  used,  however,  this  study  will  only  consider  a  flat 
bottom.    The  governing  equations  written  in  sigma  coordinates  are  as  follows: 


du  -1  dp'  1  c  dp'  3D  . 

+ +         fv       +      A_V4u 


a.    Momentum  Equations 

du  -1  dp' 

dt  p0  dx  p0  D  ca  dx 

K 

D2dG2 


dt  pn  dy  pn  D  3g  dx 


Km  <32v 

+  T3 —        + 

D2da2 


where  (7    =    z.  D 


m 


+      J\   5  +        «h(u)  (2.1) 


dv  -1  dp'  1  c  dp'  dD  , 

+ -  fu       +         A_V4v 


m 


«d(v)  (2.2) 


b.    Continuity  Equation: 

d\v  du  dv 

——      +  — —      +  —        =     o  (2.3) 

cg  dx  d\ 
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c.  Vertically  Integrated  Hydrostatic  Equation: 

p'     =     Dj^Bd^     -     f  .J  Dj^Bd^  J  d(X  (2.4) 

d.  Equation  of  State: 

B      -         ag(T-T0)  (2.5) 


e.    Thermodynamic  Equation: 

dT                                      KM  d2J 
dt                  H                   D-  co- 

_)_ 

1  £ 
— r  [S(z)-(wT')j 
D  0(5 

+ 


5d(t)  (2.6) 


In  the  above  equations,  all  horizontal  partial  derivatives  are  on  constant  sigma 
surfaces.  In  addition,  the  variables  used  are  defined  in  Table  1.  Other  symbols  in  the 
model  equations  are  listed  in  Table  2.  which  gives  values  of  constants  used  throughout 
this  study. 

2.  Domain  Size  and  Resolution 

The  model  domain  extends  6°  in  longitude  and  6°  in  latitude  (512  km  by  640 
km)  and  has  S  km  by  10  km  horizontal  resolution  with  ten  levels  in  the  vertical.  This 
is  the  finest  horizontal  resolution  that  is  practical  for  use  on  the  IBM  3033  at  the 
Naval  Postgraduate  School. 

3.  Forcing  Fields 

The  model  can  either  be  spun  up  from  rest  or  initialized  with  some  known 
current  field,  e.g..  a  baroclinic  jet  representative  of  the  mean  California  current.  In 
addition,  a  surface  heat  flux,  dependent  on  the  model  predicted  sea  surface 
temperature,  can  be  computed  and  used  to  damp  the  temperature  field  toward  a 
prescribed  climatological  value.  Details  are  given  below.  In  this  study,  the  model  is 
spun-up  using  a  representative  climatological  wind  field. 

4.  Finite  Difference  Scheme 

In  the  horizontal,  a  space-staggered  B-scheme  (Arakawa  and  Lamb.  1977)  is 
used  while  a  sigma  coordinate  system  controls  the  vertical.  The  noisefree  version  of 
the  hydrostatic  equation  in  sigma  coordinates  advocated  by  Arakawa  and  Suarez 
(1983)  has  been  implemented. 
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TABLE  1 

DEFINITIONS  OF  VARIABLES 
USED  IN  THE  MODEL 

SYMBOL 

DEFINITION 

z 

height  (positive  upwards) 

t 

time 

P' 

pressure  perturbation  from  a  vertical  average 

T 

temperature 

u,v,w 

eastward. northward,  and  vertical  (sigma) 
velocity  components,  respectively 

B 

buoyancy 

5d 

dynamic  adjustment  procedure 

S 

solar  radiation  in  the  ocean 

vv'T 

turbulent  vertical  heat  flux. 

5.  Boundary  Conditions 

The  northern,  western,  and  southern  boundaries  are  open  using  a  modified 
version  of  the  radiation  boundary  condition  of  Camerlengo  and  O'Brien  (1980).  The 
eastern  boundary,  representing  the  west  coast  of  North  America,  is  closed,  and 
straight,  with  freeor  zero-slip  capability. 

B.       SPECIFIC  EXPERIMENTAL  CONDITIONS 

1.  Method  of  Determining  the  Wind  Forcing 

The  climatological  study  of  Nelson  (1977)  was  used  to  determine  the  wind 
forcing  for  the  PE  model  domain,  which  extends  from  36.5N  to  42. 5N  and  from  124W 
to  130W  (Figure  1.1).  Nelson  has  compiled  long-term  composite  monthly  fields  of 
surface  wind  stress  on  a  1°  square  area  basis  along  the  west  coast  of  the  United  States. 
Nelson  states  that  the  monthly  mean  data  described  in  the  report  adequately  resolves 
the  seasonal  cycle,  which  is  the  dominant  time  scale  for  coastal  upwelling  (Mooers  et 
al.,  1976).  Although  the  PE  model  is  of  a  finer  horizontal  resolution  (i.e..  S  to  10  km) 
than  Nelson's  analysis  (approx.  110  km),  it  is  the  best  data  set  available  for  seasonal 
forcing  in  the  California  Current  region. 
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TABLE  2 

VALUES  OF  CONSTANTS  USED  IN  THE  MODEL 

VALUE 

NAME 

n 

2  7i  day"1 

earth  rotation  rate 

c 

p 

0.958  calgm-1  (°K)"1 

specific  heat  of  sea  water 

CD 

1.3   x   10'3 

drag  coefficient 

L 

595.  cal  gm"1 

latent  heat  of  sea  water 

T0 

278.2  °K 

constant  reference  temperature 

Pa 

1.23   x   icr3  gmcm"3 

density  of  air 

Po 

1.0276  gm  cm 

density  of  sea  water  at  TQ 

a 

2.01   x   lO"4^)"1 

thermal  expansion  coefficient 

K 

10 

number  of  levels  in  vertical 

AX 

8.   x    io5  cm 

meridional  grid  spacing 

AY 

1.   x   io6  cm 

zonal  grid  spacing 

D 

4.5   x    io5  cm 

total  ocean  depth 

<P0 

36.5°  N 

latitude  of  southern  boundary 

^m 

42.5°  N 

latitude  of  northern  boundary 

h 

124°  W 

longitude  of  eastern  boundary 

m 

130°  W 

longitude  of  western  boundary 

At 

800.  sec 

time  step 

f 

0.93   x    10'4  sec-1 

Conolis  parameter 

a 
8 

9S0.0  cm  sec"2 

acceleration  of  gravity 

AM 

2.   x    lQ17  cm4  sec"1 

biharmonic  momentum  diffusion 

AH 

2.   x  io17  cm4  sec"1 

biharmonic  heat  diffusion 

KM 

0.5  cm    sec"1 

vertical  eddy  viscosity 

Hi 

0.5  cm2  sec"1 

vertical  eddy  conductivity 

Psfc 

1013.25  mb 

surface  air  pressure 

20 


The  choice  of  selecting  a  specific  span  of  time  over  which  climatological  data 
would  be  used  was  driven  by  interest  in  the  season  most  favorable  for  coastal  upwelling 
in  the  California  Current  System,  which  has  been  most  often  observed  to  occur  during 
the  summer  (Hickey,  1979).  As  a  result,  the  months  of  June,  July,  and  August  were 
chosen  for  the  experimental  time  domain  for  this  study. 

Close  inspection  of  Nelson's  climatological  data  over  the  model's  region 
during  the  summer  months.  Figures  2.1  through  2.3.  shows  a  large  variability  in  the 
surface  stress  field  in  the  offshore  direction,  and  little  variability  in  the  alongshore 
direction.  As  a  result,  it  was  decided  to  use  an  idealized  alongshore  wind  stress  only, 
and  to  make  it  a  function  of  the  offshore  direction  only. 

A  geographical  interpretation  problem  arises  from  the  fact  that  the  PE  model 
has  a  straight  coastline.  The  climatological  data,  of  course,  is  based  on  the  actual 
position  of  the  California  coast,  which  is  not  straight.  Furthermore,  as  Figures  2.4 
through  2.6  show,  the  wind  stress  curl  has  a  strong  dependence  on  distance  from  the 
coast.  A  simple  mapping  of  data  from  the  actual  shape  of  the  coast  to  the  straight 
model  coast  was  performed  to  partly  model  the  effect  that  the  actual  coast  has  on 
values  of  wind  stress  curl. 

The  reason  for  exercising  care  in  the  preservation  of  a  realistic  pattern  of  wind 
stress  curl  is  that  wind  stress  curl  is  the  forcing  function  for  the  vertically  integrated 
mass  transport  of  the  wind-driven  ocean  circulation  (Nelson.  1977).  More  importantly, 
horizontal  variations  in  mass  transport  in  the  surface  (Ekman)  layer  produce  surface 
divergence  (convergence)  corresponding  to  positive  (negative)  wind  stress  curl.  The 
Ekman  pumping  velocity  resulting  from  the  divergence  is  approximately  (pf  )  times 
the  curl  of  the  wind  stress  (Gill,  1 9 S 2 ) .  Ekman  pumping  can  significantly  contribute  to, 
and  modify,  the  coastal  upwelling  regime. 

Another  general  observation  of  Nelson's  wind  stress  data  is  that  the  wind 
stress  vectors  during  the  months  of  June,  July,  and  August,  over  the  model  domain,  are 
parallel  to  the  coast  within  about  10-20%  accuracy.  This  means  that  the  magnitude  of 
wind  stress  vectors  represents,  to  approximately  S0-90%  accuracy,  the  alongshore 
component  of  surface  stress.  Since  the  PE  model's  coastline  is  oriented  due  north  and 
south,  the  observed  alongshore  component  is  assumed  to  coincide  with  the  model's 
north  component  of  surface  stress.  This  observation  led  to  the  decision  to  use  the 
magnitude  of  Nelson's  wind  stress  as  the  north  component  of  the  model  stress.  The 
zonal  (cross-shore)  component  of  the  model  stress  is  taken  to  be  zero  in  all 
experiments. 
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Ficure  2.1     Climatolocical  wind  stress  -  June. 
The  inside  box  outlines  the  model  domain. 
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Figure  2.2     Same  as  Figure  2.1  except  July, 
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Figure  2.4     Climatolocical  wind  stress  curl  -  June. 
"  The  inside  box  outlines  the  model  domain. 
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Figure  2.5     Same  as  Figure  2.4  except  July. 
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Figure  2.6     Same  as  Figure  2.4  except  August. 
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For  each  month  in  consideration,  a  total  of  36  data  points  were  extracted 
from  Nelson's  data  set  for  both  wind  stress  and  wind  stress  curl  because  the  domain  of 
the  PE  model  spans  6°  of  latitude  and  6°  of  longitude  and  Nelson's  data  is  of  a  1° 
resolution.  Since  the  alongshore  variation  was  considered  much  less  significant  than 
the  cross-shore  variation  for  both  wind  stress  and  wind  stress  curl,  an  alongshore 
average  of  the  wind  stress  curl  was  calculated  from  Figures  2.4  through  2.6  and  then 
averaged  over  the  three  months  of  June,  July,  and  August.  The  result  is  shown  by  the 
circles  connected  by  a  dotted  line  in  Figure  2.7.  The  dominant  feature  is  a  region  of 
strong  positive  curl  within  200  km  of  the  eastern  boundary,  and  weaker  negative  curl 
over  the  outer  300  km.  This  time-averaged  and  alongshore  averaged  curl  reaches  its 
maximum  value  at  the  eastern  boundary. 

The  modeling  effort  for  producing  an  idealized  wind  forcing  function  was 
based  on  understanding  that  the  wind  stress  curl  is  a  prominent  feature  which  will 
produce  a  modifying  effect  on  the  dynamics  of  ocean  circulation.  The  observed  three- 
month  average,  alongshore-averaged  wind  stress  curl,  derived  from  Nelson's 
climatology  (Figure  2.7),  was  closely  examined  for  features  which  could  be  represented 
by  an  analytical  function.  A  Gaussian  function  in  the  offshore  coordinate  was  found 
to  fit  the  climatological  data  satisfactorily.  The  Gaussian  function  used  to  represent 
the  observed  curl  is: 

dx/dx  =  0.77exp[-0.5(XR/95.0)2)  -.11  ,  (2.7) 

where  XR  is  distance  offshore  in  kilometers.  Figure  2.7  shows  a  plot  of  this  analytical 
function  (square  boxes  connected  by  a  solid  line)  along  with  the  observed  values,  and 
the  fit  is  seen  to  be  quite  good. 

The  PE  model  is  designed  to  use  wind  input  (for  convenience  of  wind  stress 
and  heat  flux  computations)  rather  than  wind  stress  curl,  but  in  order  to  benefit  from  a 
realistic  wind  stress  curl  function  and  also  retain  the  ability  to  force  the  PE  Model  with 
wind  data,  idealized  wind  data  were  derived  from  the  above  Gaussian  representation  of 
wind  stress  curl  in  the  following  manner.  The  first  step  was  to  calculate  the  value  of 
the  wind  stress  curl  (at  each  u,v  gridpoint)  using  the  Gaussian  function  in  Figure  2.7. 
Next,  to  compute  the  surface  stress  (at  the  temperature  gridpoints)  from  the  wind 
stress  curl  (at  each  u,v  gridpoint),  the  wind  stress  curl  was  numerically  integrated  in  the 
offshore  direction.  The  constant  of  integration  was  obtained  from  the  alongshore  and 
three-month  average  of  the  data  from  Nelson  (1977).     This  constant  is   t     =  -0.65 
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dyne; cm"  and  it  represents  the  value  of  the  alongshore  wind  stress  at  the  coast.  The 
final  step  in  deriving  winds  from  the  idealized  stress  was  to  utilize  the  bulk 
aerodynamic  formula  to  compute  winds  at  model  temperature  gridpomts.  A  constant 
drag  coefficient  of  0.0013  was  utilized  with  an  air  density  of  0.00123  gmcnr.  Wind 
values  were  defined  at  model  temperature  gridpoints  because  the  surface  wind  speed  is 
used  in  heat  budget  calculations  within  the  PE  model,  and  therefore  it  is  needed  at  each 
temperature  gridpoint.  Upon  completion  of  heat  budget  calculations,  the  PE  model 
computes  the  wind  stress  and  performs  a  four-way  averaging  to  get  the  stress  at  each 
u.v  gridpoint,  whereupon  it  drives  model  surface  currents. 

A  plot  of  the  wind  stress,  obtained  by  integration  of  the  idealized  curl 
function,  is  shown  by  the  squares  in  Figure  2.S.  It  is  clear  that  the  proper  shape  and 
order  of  magnitude  has  been  captured  well.  Any  difference  between  modeled  and 
climatological  wind  stress  is  probably  due  to  the  modeling  assumption  that  the  curl  is 
determined  solely  by  the  offshore  variation  of  the  alongshore  stress.  A  plot  of  the 
alongshore  component  of  the  winds  computed  from  the  observed  and  modeled  stress, 
respectively,  using  the  bulk  aerodynamic  formula,  is  shown  in  Figure  2.9.  Again,  the 
shape  and  magnitude  is  well  represented  and  any  differences  are  just  reflections  of 
differences  in  the  wind  stress  plots.  This  idealized  wind  (having  the  stress  and  stress 
curl  as  shown  in  Figures  2.8  and  2.7  respectively)  is  the  actual  data  which  was  used  to 
force  the  PE  model  in  some  of  the  wind  forcing  experiments. 
2.   PE  Model  Surface  Thermal  Damping 

Surface  thermal  damping  becomes  a  crucial  requirement  for  realistic  modeling 
of  ocean  temperature  fields  in  the  presence  of  atmospheric  surface  wind  forcing.  This 
is  because  the  summer  heat  flux  in  the  California  region  is  a  downward  flux  due  to  a 
greatly  reduced  sensible  and  latent  heat  flux  produced  by  the  existence  of  cold 
(upwelled)  water.  The  stronger  the  alongshore  (southward)  winds,  the  greater  the 
upwelling  and  the  colder  the  surface  water  and  the  greater  the  downward  heat  flux. 
The  downward  heat  flux  is  thus  negatively  correlated  with  the  sea  surface  temperature, 
hence  it  is  a  damping  factor.  Another  way  to  understand  the  damping  (as  contrasted 
to  forcing)  nature  of  the  surface  heat  fluxes  in  the  CCS  is  to  recognize  that  the  largest 
downward  surface  heat  fluxes  occur  where  the  ocean  surface  temperatures  are  the 
coldest.  Heat  flux  variations  related  to  atmospheric  variations  in  clouds  or  air 
temperature  are  not  considered  to  be  an  important  forcing  factor  in  the  CCS.  Heat 
fluxes    are    therefore   computed    from   bulk    formulae   using    prescribed   quantities    of 
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atmospheric  solar  radiation,  clouds,  surface  air  temperature,  relative  humidity  and 
winds  (Haney  et  al,  1978).  The  following  representative  mean  values,  for  the  summer 
season,  are  used:  fractional  cloud  coverage  of  48%,  solar  insolation  at  top  of  the 
atmosphere  of  983.3  cal  cm  /day,  and  surface  relative  humidity  of  85%.  The  value  of 
solar  radiation  was  obtained  from  the  Smithsonian  Meteorological  tables  (List,  1963), 
the  fractional  cloud  cover  was  obtained  from  the  satellite  cloud  atlas  of  Miller  and 
Feddes  (1971),  and  the  relative  humidity  was  computed  from  climatological  data  of 
Nelson  (19S3). 

Sensible  and  latent  heat  fluxes  are  calculated  using  exchange  coefficients  which 
are  proportional  to  the  wind-speed-dependent  neutral  drag  coefficient  [C^(V)j  of  Large 
and  Pond  (1981).  Initial  conditions  for  sea  surface  and  air  temperatures,  respectively, 
are  15°  C  sea  surface  temperature  and  13°  C  air  temperature.  The  sea  surface 
temperature  used  is  a  representative  offshore  temperature  for  the  summer,  whereas  the 
air  temperature  was  chosen  so  that  there  would  initially  be  no  net  heat  flux  at  the 
surface  with  a  steady  wind  speed  of  ~  850  cm  sec  (wind  speed  used  for  experiment  1). 
The  imposed  criteria  of  no  initial  heat  flux  guaranteed  that  any  heat  flux  generated 
after  model  initialization  was  only  in  response  to  the  changing  sea  surface  temperature. 
This  criteria  also  assured  that  the  model'  heat  budget  acted  only  to  thermally  damp  the 
effects  of  surface  wind  forcing  and  not  as  an  additional  forcing  function.  This  is  not 
only  realistic,  as  noted  above,  but  it  also  simplifies  the  interpretation  later. 

The  damping  nature  of  the  computed  heat  flux  is  demonstrated  best  by 
considering  a  linearized  form  of  the  modeled  equations  for  the  total  downward  surface 
heat  flux.  Q.  which  is  Q  =  y(TA  -SST).  TA  is  an  apparent  atmospheric  equilibrium 
temperature,  SST  the  model  surface-layer  temperature  and  y  a  coupling  coefficient 
strongly  dependent  on  surface  wind  speed,  (Haney,  1985).  The  corresponding 
relaxation  time  for  the  surface  temperature  is  given  by  X"1  =  pC  hy"1  (Haney.  1985), 

r 

where  p  is  the  density  of  sea  water.  C  the  specific  heat  of  sea  water,  and  h  the  depth 
of  the  mixed  layer.  With  mean  climatological  winds.  y~25Wm"2K"1  (Haney,  1971). 
Using  h  =  50m,  a  value  is  obtained  for  \~l  of  100  days.  Since  y  is  directly 
proportional  to  wind  speed,  if  the  wind  speed  changes  on  the  order  of  20%,  the 
temperature  relaxation  time  will  also  change  by  20%  yielding  a  range  for  X"1  of  SO  to 
120  days.  This  long  damping  time  will  not  mask  any  change  in  ocean  temperature  due 
to  dynamical  processes  which  are  expected  to  have  a  time  scale  of  only  several  days. 
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III.  RESULTS  OF  WIND  EXPERIMENTS 

A.       EXPERIMENTS  WITH  FREE  SLIP  BOUNDARY  CONDITIONS 
1.  Experiment  1  (Uniform  Wind) 

In  the  first  experiment  conducted,  the  idealized  winds  obtained  from  the 
procedure  described  in  Chapter  II.  and  shown  in  Figure  2.9,  were  averaged  in  the 
cross-shore  direction  resulting  in  a  constant  value  (approximately  850  cm  sec).  The 
forcing  for  experiment  1  was,  therefore,  uniform  (alongshore  and  cross-shore  invariant) 
and  steady  (time  invariant).  This  idealized  wind  was  considered  to  be  the  least 
complex;  yet  it  still  retains  some  of  the  realism  of  observed,  long-term,  mean 
climatological  winds  during  the  summer  season  over  the  geographical  region 
encompassed  by  the  PE  model  domain.  All  experiments  in  this  study,  unless  otherwise 
stated,  included  surface  thermal  damping  as  described  in  Chapter  II.  Additionally,  the 
eastern  boundary  (representing  the  west  coast  of  North  America)  was  set  as  a  free-slip 
condition  in  this  experiment,  meaning  that  the  tangential  ocean  velocities  were 
unrestricted  by  frictional  effects. 

At  this  time,  it  is  important  to  discuss  the  initial  conditions,  in  particular  the 
model  temperature  stratification  for  the  wind  forcing  experiments.  A  pure  exponential 
temperature  profile,  as  a  function  of  height  (z),  was  used  to  initialize  all  of  the 
experiments  within  this  study.  This  profile  has  a  length-scale  of  450  m  with  a 
temperature  at  the  first  sigma  level  (ocean  surface  level)  of  15°  C  and  a  value  at  great 
depth  of  2°  C.   The  exact  form  of  this  temperature  function  is  as  follows: 

T(z)    =    2+13  ez45°.  (3.1) 

The  temperature  profile  was  derived  from  observations  used  to  support  the 
Princeton  Dynalysis  model  (Blumberg  er  al.,  1984)  and  is  representative  of  the  long- 
term,  mean  climatological  temperature  stratification  for  the  California  coastal  region. 
Additionally,  the  Brunc-Vaisala  frequency  profile,  N  (z),  was  calculated  analytically 
using  the  temperature  function  in  (3.1)  from 

N2    -     a  2  <?T  <3z  .  (3.2) 
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The  model  temperature  profile  and  resultant  N  profile  are  shown  in  Figures  3.1  and 
3.2.  Since  the  initial  temperature  field  was  horizontally  uniform,  the  initial  ocean 
currents  were  set  to  zero. 

The  model  was  spun  up  by  turning  on  the  wind,  at  model  time  zero,  to  its  full 
magnitude.  The  impulsive  wind  produced  an  impulsive  surface  stress  which  resulted  in 
inertial  oscillations  of  the  near-surface  ocean  currents.  These  oscillations  are 
demonstrated  by  a  time  series  of  the  v-component  of  velocity  at  a  selected  gridpoint 
away  from  the  coast  (Figure  3.3).  The  inertial  period  approximated  from  Figure  3.3  of 
~  0.75  day  compares  well  with  the  theoretical  inertial  period  of  2tt  f.  Off-shore 
Ekman  transport  was  apparent  from  the  negative  values  of  the  u-component  of 
velocity  throughout  the  model  spin-up  (see  Figure  3.11).  Since  the  wind  was  northerly, 
and  the  model  domain  is  in  the  northern  hemisphere,  Ekman  transport  was  directed  to 
the  right  of  the  direction  of  the  surface  stress;  which  is  towards  the  west  (represented 
by  negative  values  of  the  u-component  of  velocity).  The  Ekman  flow  away  from  the 
rigid  eastern  boundary  produces  some  very  important  effects.  Water  carried  seaward, 
in  the  oceanic  boundary  layer,  produces  divergence  at  the  coast.  Due  to  the  principle 
of  conservation  of  mass,  the  void  created  by  coastal  divergence  must  be  filled  by 
upwelling  of  colder  water  from  lower  levels. 

The  effect  of  upwelling.  which  brings  colder  water  up  from  the  lower  levels,  is 
to  perturb  the  initial,  horizontally  uniform,  temperature  structure.  The  isotherms  are 
bent  upward  near  the  coast,  which  establishes  a  horizontal  temperature  gradient  m  the 
upper  levels  of  the  water  column.  The  vertical  velocities  established  by  upwelling  in 
the  PE  model  compare  favorably  with  an  analytical  approximation  of  upwelling  rates 
using  basic  principles  and  assumptions  (see  Appendix).  Consistent  with  geostrophy  (a 
balance  between  the  pressure  gradient  force  and  the  Coriolis  force)  and  the  resulting 
thermal  wind  equations  (a  relation  of  vertical  shear  to  the  horizontal  temperature 
gradient),  the  model  spin-up  is  characterized  by  a  coastally-trapped,  alongshore  jet  in 
the  upper  1000  m  of  the  water  column.  A  simple  schematic  of  the  coastal  jet  is  shown 
in  Figure  3.4  from  Gill  (19S2).  The  southward  flowing  jet  along  the  eastern  boundary 
of  the  model  domain  will  be  the  focal  point  of  this  study. 

Figures  3.5  through  3.7  show  contours  of  the  v-component  of  surface  velocity, 
along  with  contours  of  temperature,  during  the  spin-up  phase  at  days  30,  60,  and  90 
respectively.  The  coastal  jet  develops  in  response  to  the  increase  in  the  horizontal 
temperature    gradient    with    time.     The    near-coast,    surface,    horizontal    temperature 
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Figure  3.1     Model  initial  temperature  profile. 
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Figure  3.3     V-component  inertial  oscillations. 
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Figure  3.4     Schematic  of  coastal  jet  formation  (from  Gill,  19S2). 
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gradient  at  day  30  (Figure  3.5)  is  0.055°/km  corresponding  to  a  maximum  surface, 
southward  alongshore  velocity  component  of  60  cm  sec.  At  day  60  (Figure  3.6),  the 
gradient  of  temperature  has  increased  to  0.063°  km  with  an  associated  maximum 
southward  velocity  of  SO  cm  sec.  Finally,  at  day  90  (Figure  3.7),  the  gradient  of 
temperature  has  further  increased  to  0.070°/km  with  a  maximum  southward  velocity  of 
100  cm;  sec. 

The  vertical  structure  of  the  developed  jet  is  best  described  by  computing  the 
alongshore  and  ten-day  time-average  of  the  v-component  of  velocity  and  plotting  that 
field  as  a  function  of  cross-shore  direction  and  depth.  Figure  3.S  shows  some 
interesting  features.  The  depth  of  penetration  of  the  southward  flowing,  coastal  jet  is 
approximately  1000  m  with  an  offshore  extent  of  approximately  100  km.  The 
maximum  surface  velocity  lies  approximately  25  km  offshore  and  the  jet  axis  is 
vertically  oriented  with  depth.  The  first  internal  Rossby  radius  of  deformation  was 
calculated  to  be  29.3  km.  by  the  method  of  Feliks  (19S5).  which  coincides  well  with  the 
position  o(  the  surface  maximum,  alongshore  velocity  component  in  this  experiment. 
An  undercurrent  is  present  which  is  characterized  as  a  deep,  slow  (~  5  cm  sec), 
poleward-flowing  current  underlying  the  southward  flowing  jet.  This  undercurrent  has 
a  maximum  alongshore  velocity  component  at  the  coast.  Additionally,  note  that  the 
weakest  vertical  shear  occurs  on  the  anticyclonic  side  (offshore  or  right  side)  of  the 
southward-flowing  jet.  While  it  may  not  be  easy  to  see  from  the  figure,  the  largest 
magnitude  of  relative  vorticity  also  occurs  on  the  anticyclonic  side  of  the  jet  in  this 
experiment. 

A  cross-section  of  the  alongshoreand  time-averaged  temperature  field  is  shown 
in  Figure  3.9.  Away  from  the  coast,  the  initial  temperature  stratification  o[  15"  C  at 
the  surface  and  2°  C  at  depth,  has  only  been  slightly  modified,  if  at  all.  However,  the 
coastal  stratification  demonstrates  a  classical  upwelling  regime  by  the  upward  bending 
isotherms.  As  previously  described,  the  horizontal  temperature  gradients  established 
support  the  geostrophic  jet  (Figure  3.8).  The  surface,  coastal  temperature  has  been 
cooled,  from  an  initial  temperature  of  15°  C,  to  approximately  9°  C.  This  coastal 
temperature  structure  has  important  implications  in  potential  vorticity  analysis  and  will 
be  discussed  later. 

2.  Experiment  2  (Wind  with  Curl) 

Experiment  2  differed  from  experiment  1  only  in  that  the  forcing  function 
included  the  important  effect  of  wind  stress  curl.     Idealized  wind  data  (Figure  2.9) 
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Figure  3.7     Same  as  Figure  3.5  except  day  90. 


43 


WW  0\  :■ 


o 
-o 


oooocooooooooooooooooooooooooooooooooooooooooo 

CsirNJC\IC\ICMCNOJC\|(NC\ICNC\IOJCNCNC\IC\ICsl(N0^CNC\l<NC\ICNC\;CNrM 

1    '        M    M    '    I    '    M    I    I    I    I    I    I    |    I    |    |    I    |    |    |    |    |    |    I    |    |    |    |   |    |    |    |    |    I    |    |    |    I 


(uu)  M^dsQ 


Figure  3.S     Cross-section  of  v-component  velocitv 
walonsshore-  and  time  averaged  (aavs  8  1-90) 
expenment  1  -  contour  interval  =  3  cm  sec. 


44 


i  ! 


i  CM 


CO 


>o 


o 
■o 
in 


o 
■m 


o 
■o 


o 

■in 
m 


c 

o  .^ 

o  "— 

^     CD 

u 

c 
o 


o 
n 


o 


^n 


en 

'"D 
O 

CD 


CO 
SO 


oooooooocoooooooooooooooooooooooocoooooooooooo 

CM  CM  CM  CM  CM  CM  CM  C  J  O  J  CM  CM  CM  CM  CM  CM  CM  C\!  fN  CM  CM  CM  CM  CM  CM  CM  CM  O  J  CM  CM  CM  IN  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM 
|V^^^'n^^CCCT)Q^rJ'Ov-ini0^cOC7>0^1CNrO-^uO(Or^cOC7)Q^CMrO'^ir.cOI^CDaiO^CMfO-*iO 

'    I    M    I    I    |    I    |    |    I    |    |    I    |    |    |    I    |    |    |    |    |    |    |    |    |    I    |    |    |    I 


(iu)  q^dsa 


Figure  3.9     Same  as  Figure  3.S  except  temperature 


contour  interval  =   1 


E 


45 


derived  from  the  Gaussian  representation  of  wind  stress  curl  (eqn.  2.7),  as  described  in 
Chapter  II.  was  used  to  force  the  PE  model  in  experiment  2.  The  inclusion  of  wind 
stress  curl  advanced  the  complexity,  and  realism,  of  model  forcing  because  the  wind 
was  no  longer  uniform,  as  was  the  case  in  experiment  1.  Otherwise,  all  aspects  of  the 
first  two  experiments  were  identical. 

The  idealized  wind  stress  curl  (Figure  2.7)  has  a  maximum  value  at  the  coast 
of  0.65  dynes/ cm"/ 100  km  which  decreases  to  0.00  dynes; cm-' 100  km  at  a  distance  of 
187  km  offshore.  After  the  zero  crossing,  the  wind  stress  curl  continues  to  slowly 
decrease  to  a  value  of  -0.11  dynes  cm2 TOO  km  at  the  western  boundary  of  the  PE 
model  domain. 

Nelson  (1977)  described  the  observed,  long-term,  mean  climatological  wind 
stress  curl  off  the  California  coast  as  follows: 


Existence  of  an  offshore  wind  stress  maximum  results  in  a  line  of  zero  wind  stress 
curl  approximately  parallel  to  the  coast.  Positive  curl  occurs  inshore  of  the 
maximum  wind  stress.  Negative  curl  in  the  offshore  resion  is  associated  with  the 
anticyclonic  atmospheric  circulation  over  the  interior  "ocean.  The  positive  curl 
near  "the  coast  is  related  to  toposraphv  and  to  local  features  in  the  surface  wind 
stress  distributions. 


The  main  effect  caused  by  wind  stress  curl  is  that  the  variation  of  wind  stress 
in  the  offshore  direction  causes  a  variation  in  the  Ekman  transport  which  leads  to 
convergence  and  divergence  of  mass.  With  the  northerly  wind  blowing  along  the 
eastern  boundary  of  the  model,  the  positive  curl  produces  oceanic,  upper-level, 
horizontal  divergence  so  fluid  from  below  must  be  "sucked"  vertically  into  the 
boundary  layer  to  replace  fluid  transported  away  by  mass  divergence.  The  effect  is 
Ekman  pumping.  Negative  curl,  conversely,  creates  mass  convergence  and  causes 
downwelling  (Figure  3.10).  As  previously  discussed,  the  wind  stress  curl  goes  through 
zero  at  187  km  offshore  and,  therefore,  the  positive  wind  stress  curl  along  the  coast, 
causing  Ekman  pumping,  will  have  a  modifying  effect  on  the  coastal  upwellmg  process 
described  in  section  III.A.l.  The  effect  of  an  offshore  varying  wind  stress  curl,  as  in 
experiment  2,  is  depicted  by  contours  of  the  u-component  of  velocity  during  the  initial 
model  spin-up  phase.  Figure  3.11  clearly  demonstrates  that  Ekman  transport  is  a 
function  of  the  strength  of  the  northerly  winds.  Referring  back  to  Figure  2.9,  it  is 
noted  that  the  wind  increases  offshore  to  approximately  190  km  and  Figure  3.11 
displays  the  corresponding  increase  in  magnitude  of  the  u-component  to  approximately 
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190  km.  The  u-component  is  negative  (offshore)  everywhere  within  the  domain  which 
is  also  consistent  with  Ekman  transport  being  directed  to  the  right  of  the  wind.  This 
offshore  increase  in  the  magnitude  of  the  u-component  of  velocity  is  a  representation 
of  mass  divergence  resulting  in  Ekman  pumping.  Accordingly,  proceeding  past  the 
zero  crossing  of  wind  stress  curl  further  offshore,  the  magnitude  of  the  u-component 
contours  begins  to  decrease  which  is  an  indication  of  mass  convergence. 

A  cross-section  of  the  alongshore-  and  time-averaged  v-component  of  velocity 
from  experiment  2  is  displayed  in  Figure  3.12.  Striking  differences  in  the  structure  and 
magnitude  of  the  jet  are  immediately  discernible  in  this  figure  when  compared  with  a 
similar  plot  from  experiment  1  (Figure  3.S).  The  offshore  extent  of  the  surface  coastal 
jet  in  experiment  2  has  increased  by  almost  90  km  as  compared  with  experiment  1. 
Note,  however,  that  the  effect  is  very  shallow  and  quite  weak.  Albeit  its  relatively 
weak  effect,  the  vertical  shear  on  the  offshore  (anticyclonic)  side  of  the  coastal  jet  is, 
indeed,  strengthened  over  that  in  experiment  1.  The  most  noticable  difference  between 
the  two  profiles  is  the  difference  in  the  strengths  of  the  jet.  Figure  3.S  shows  the 
maximum  surface  velocity  of  the  jet  to  be  -105  cm  sec  whereas  Figure  3.12  shows  the 
maximum  to  be  only  -75  cm  sec.  The  explanation  ot~  this  lies  with  the  value  of  the 
wind  at  the  coast  in  each  experiment.  Experiment  1  used  a  constant  value  of 
approximately  -S50  cm/sec  (domain  climatological  value)  whereas  the  coastal  value  of 
wind  in  experiment  2  was  only  about  -630  cm  sec  (coastal  climatological  average 
value).  Coastal  upwelling  induced  by  the  wind  stress  near  the  coast  is  clearly  the 
dominant  factor  over  Ekman  pumping.  A  difference  of  only  2  m  sec  (  ~  5  kts)  average 
wind  speed  can  obviously  produce  a  dramatic  change  in  the  oceanic  response  adjacent 
to  the  coast  of  California  when  the  wind  is  northerly. 


B.       EXPERIMENTS  WITH  ZERO  SLIP  BOUNDARY  CONDITIONS 

1.  Experiment  3  (Uniform  Wind) 

The  zero-slip  boundary  condition  at  the  eastern  boundary  was  introduced  into 
experiments  3  and  4.  The  importance  of  the  free-slip  versus  the  zero-slip  boundary 
condition  was  addressed  in  Chapter  II  and  it  has  been  the  subject  of  study  by  ocean 
circulation  modelers  since  the  studies  of  Bryan  (1963)  and  Veronis  (1966).  Blandford 
(1971)  demonstrated  that  the  large  difference  between  solutions  obtained  by  Veronis 
and  Bryan  was  caused  by  the  difference  between  the  free-slip  and  zero-slip  conditions 
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Figure  3.10     Schematic  of  Ekman  pumping  (from  GUI,  1977). 
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Fisure  3.11     Contours  of  u-component  velocity 

contour  interval  =  0.5  cm  sec 

experiment  2  -  day  6. 
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Figure  3.12     As  in  Figure  3.8  except  experiment  2. 
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used.  Veronis'  free-slip  condition  produced  a  maximum  velocity  at  the  coast  while  the 
zero-slip  condition  used  by  Bryan  produced  a  maximum  velocity  offshore.  It  is  known 
that  barotropic  instability  can  occur  if  the  profile  curvature  changes  sign  (Haltiner  and 
Williams,  19S0)  and  the  zero-slip  solution  of  Bryan's  certainly  satisfied  that  necessary 
condition. 

Experiment  3,  as  stated  previously,  included  the  zero-slip  boundary  condition. 
All  other  aspects  of  the  model  were  identical  to  experiment  1.  Figure  3.13  is  a 
normalized  plot  of  the  alongshore  (v)  component  of  velocity  averaged  in  time  (days 
S6-90).  in  depth  (surface  to  ~  1000  m),  and  alongshore.  For  comparison,  a  similar 
plot  is  displayed  for  experiment  1  (Figure  3.14).  The  effects  of  the  zero-slip  boundary 
condition  employed  in  experiment  3  are  immediately  apparent.  The  frictional  influence 
of  zero-slip  on  the  near-shore  tangential  velocity  (v-component)  dramatically  reduced 
the  coastal  velocity  as  compared  to  the  coastal  velocity  generated  in  experiment  1  with 
the  free-slip  condition.  Therefore,  the  main  effect  of  the  zero-slip  boundary  condition 
is  to  significantly  increase  the  cyclonic  horizontal  shear  on  the  inshore  side  of  the  jet. 
The  cyclonic  horizontal  shear  in  experiment  1  was  .15  x  10  s"1  whereas  the  zero-slip 
condition  increased  the  shear  in  experiment  3  to  .6  x  10  s  .  Another  effect  of  the 
zero-slip  condition  in  experiment  3  was  to  move  the  depth-averaged  coastal  jet  axis 
approximately  8  km  farther  offshore  than  that  in  experiment  1  (free-slip). 

As  discussed  earlier,  a  major  result  of  Bryan's  zero-slip  boundary  condition 
study  was  that  the  maximum  alongshore  velocity  was  located  offshore  causing  the 
necessary  condition  of  barotropic  instability  (profile  curvature  changes  sign)  to  be 
satisfied,  whereas  the  free-slip  study  by  Veronis  did  not  satisfy  that  necessary  condition. 
In  this  study,  the  zero-slip  condition  was  not  an  essential  requirement  to  satisfy  the 
necessary  condition  of  barotropic  instability.  It  can  be  seen  that  experiment  1,  which 
had  the  free-slip  condition,  also  satisfied  the  necessary  condition  (evident  by  the  jet 
profile  curvature  changing  sign  approximately  75  km  offshore  in  Figure  3.14).  The 
upwellmginduced  coastal  jet,  simulated  in  this  eastern  boundary  study,  has  a  velocity 
maximum  offshore  regardless  of  the  boundary  condition  used  (free-slip  or  zero-slip). 
In  contrast,  the  western  boundary  current  in  the  free-slip  study  by  Veronis  had  a 
maximum  velocity  at  the  coast  which  prevented  the  profile  curvature  from  changing 
signs.  Apparently,  the  baroclinicity  of  the  modeled  upwelling-induced  coastal  jet 
dominates  the  positioning  of  the  jet  axis  over  the  effects  of  the  frictional  eastern 
boundary  condition. 
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Figure  3.14     Same  as  Figure  3.13  except  experiment  1, 
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It  is  again  interesting  to  analyze  the  timeand  alongshore-averaged  cross- 
section  of  the  v-component  of  velocity.  Figure  3.15  is  the  profile  generated  by 
experiment  3.  The  effect  of  a  zero-slip  boundary  condition  is  immediately  apparent  in 
this  plot.  The  tangential  velocities  along  the  coast  in  experiment  3  are  much  less  than 
those  generated  in  experiment  1  (Figure  3.8).  Because  of  this,  the  coastal  isotachs.  in 
experiment  3  are  sloped  upward  which  significantly  reduces  the  vertical  shear  on  the 
inshore  (cyclonic)  side  of  the  jet  as  compared  with  experiment  1.  The  result  is  that  the 
jet  is  more  symmetric  in  the  offshore  direction.  However,  the  positive  voracity  on  the 
inshore  side  of  the  jet  is  new  about  twice  as  big  as  the  negative  vorticity  on  the 
offshore  side.  This  is  very  different  from  the  pattern  in  experiment  1  (Figure  3.8). 
Another  interesting  result  of  the  zero-slip  boundary  condition  is  that  the  northward 
flowing  undercurrent  maximum  has  also  moved  further  offshore  along  with  the  coastal 
jet. 

2.  Experiment  4  (Wind  nith  Curl) 

As  an  extension  of  the  zero-slip  boundary  condition  study,  experiment  4 
incorporated  zero-slip  on  the  eastern  boundary  along  with  the  wind  stress  curl  as  in 
experiment  2.  Except  for  the  zero-slip  condition,  all  other  aspects  of  the  model  in 
experiment  4  were  identical  to  experiment  2.  As  would  be  expected,  the  basic 
characteristics  of  the  zero-slip  boundary  condition,  as  described  in  the  results  of 
experiment  3,  were  reproduced  in  experiment  4.  Again,  the  vertical  shear  on  the 
cyclonic  side  of  the  jet  was  decreased  by  the  upward-bending  isotachs  near  the  coast. 
As  in  experiment  3,  this  increased  the  horizontal  shear,  and  hence  the  relative  vorticity 
on  the  inshore  side  of  the  jet.  The  undercurrent  maximum  is  the  narrowest  of  all  four 
experiments  and  its  alongshore  velocity  maximum  lies  furthest  offshore  due  to  the 
combined  effects  of  wind  stress  curl  and  zero-slip.  Figure  3.16  demonstrates  the  salient 
features  of  the  mean  velocity  profile. 

C.       STABILITY  ANALYSIS 

It  is  now  well  known  that,  at  times,  the  instantaneous  California  Current  consists 
of  "intense  meandering  current  filaments  (jets)  intermingled  with  synoptic-mesoscale 
eddies"  (Mooers  and  Robinson,  1984).  The  exact  generation  mechanism,  or  source,  of 
the  observed  eddies  and  their  role  in  the  dynamics  of  the  CCS  have  yet  to  be 
determined.  One  hypothesis  currently  under  investigation  is  that  the  eddies  are  locally 
generated  through  dynamic  instability  of  the  seasonal  mean  flow.    Thomson  (1984) 
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Figure  3.15     Same  as  Figure  3.8  except  experiment  3. 
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Figure  3.16     Same  as  Figure  3.S  except  experiment  4. 
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presented  a  strong  case  of  an  observed  cyclonic  eddy  which  formed  over  the 
continental  margin  of  Vancouver  Island  through  the  process  of  baroclinic  instability 
with  an  additional,  yet  secondary,  contribution  from  barotropic  instability.  He 
furthermore  suggested  that  the  energy  for  the  instability  was  derived  from  the  vertical 
shear  of  the  geostrophic  flow  associated  with  a  southeastward  surface  current  and  a 
northwestward  subsurface  flow  representative  of  the  California  Undercurrent. 

In  the  previous  sections  of  this  chapter,  the  coastal  jet  profile  was  described  for 
each  experiment  and  was  shown  to  satisfy  the  necessary  condition  for  barotropic 
instability.  The  results  of  Thomson  (1984)  suggest  that  a  more  important  source  may 
be  baroclinic  instability  of  the  mean  flow.  Of  the  four  experiments  discussed  in  this 
chapter,  none  developed  instability  and,  therefore,  eddies  and  current  filaments  were 
not  generated.  However,  it  is  extremely  important  to  analyze  the  mean  How  to 
establish  its  potential  for  mixed  (barotropic  and  baroclinic)  instability.  If  potential  for 
instability  indeed  exists  in  the  mean  How  generated  by  the  PE  model,  it  may  provide 
insight  into  the  complexities  of  eddy  generation  in  the  CCS. 

Implications  of  the  distribution  of  potential  vorticity  as  a  signature  of  instability 
have  been  examined  by  Watts  and  others  (Robinson,  1983)  for  the  Gulf  Stream.  Watts 
stated  that  potential  vorticity,  used  as  a  conservative  quantity,  can  serve  as  a  tracer  in 
the  ocean.  Additionally,  the  cross-stream  distribution  of  potential  vorticity  may  be 
used  as  a  test  of  predictions  by  quasigeostrophic  theory  regarding  necessary  conditions 
for  unstable  growth  of  fluctuations.  The  conditions  which  are  necessary  for  instability 
to  occur  are  very  useful  because,  if  they  are  not  satisfied,  it  can  be  concluded  that  the 
dynamical  constraints  will  not  allow  available  energy  to  be  released  (Gill,  1982).  One 
of  the  conditions  is  that  the  cross-stream  derivative  of  potential  vorticity  must  change 
sign  somewhere  within  the  domain.  Another  condition  is  that  the  product  of  the  cross- 
stream  derivative  with  the  basic  current  must  be  positive.  Also,  theory  formally 
requires  a  basic  state  current  which  is  slowly  changing  m  space  and  time,  so  that 
linearization  applies  for  perturbation  quantities  (Robinson,  1983)  ;  this  requirement  is 
met  by  the  structure  of  the  coastal  jets  generated  in  this  study. 

The  following  expression,  in  cartesian  coordinates,  used  by  Watts  to  examine  the 
potential  vorticity  (q)  signature  in  the  Gulf  Stream,  will  by  used  to  examine  potential 
instability  of  the  coastal  jet  generated  in  this  study: 
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Figure  3.17  shows  a  cross-section  of  contours  of  the  time-averaged  (days  81-90) 

and  meridionally-averaged  potential  vorticity  for  experiment  1  in  units  of  :C  m     sec"1 

scaled  by  106.    Similarities  in  the  structure  of  the  potential  vorticity  field  between  the 

Gulf  Stream  (Figure   3.18)  and  the   simulated  California  Current  (Figure   3.17)  are 

apparent.    There  is  a  tendency  for  potential  vorticity  to  be  uniform  along  isothermal 

surfaces  and  also  to  change  vertically  in  the  same  manner  as  the  offshore  temperature 

stratification.    The  potential  vorticity  field  for  experiment  1  is  further  characterized  by 

the  following  features: 

The  range  of  potential  vorticity  lies  between  0.0-2.5   x   10"6  DC  rri     sec"1. 

A  relative  minimum  exists  in  the  surface  laver  caused  bv  weak  stratification  due 
to  turbulent  vertical  mixing. 

A  stronser  relative  minimum  exists  in  a  thicker  surface  laver  near-shore  due  to 
upwelling  induced  weak  stratification. 

A  relative  maximum  in  potential  vorticity  exists  in  the  middle  to  western  part  of 
the  domain  at  a  depth  or  ~  90  m  (the  "seasonal"  thermocline  in  the  model). 

A  change  in  the  sign  of  the  slope  of  potential  vorticity  contours  is  apparent  in 
the  east-west  direction  of  the  upper  levels,  near-shore. 

One  difference  between  the  potential  vorticity  field  in  the  Gulf  Stream  and  that  in 

the  simulated  CCS  is  the  difference  in  magnitude  of  potential  vorticity.    The  magnitude 

of  potential  vorticity  in  the  Gulf  Stream  is  approximately  twice  that  in  the  simulated 

CCS.     This    difference    is    mainly    attributable    to    the    difference    in    temperature 

stratification  between  the  Gulf  Stream  and  the  CCS.    From  the  isotherms  depicted  in 

Figure    3.18    it    can   be    seen    that    the    strongest    stratification   in    the    Gulf  Stream 

approaches  "~   10  °  C  over  140  m  or  0.07  °C  m.  whereas  from  Figure  3.9  the  maximum 

temperature  stratification  for  experiment  1  is  only  ~  3.5  °C  over  1-40  m  or  0.03  DC  m. 

Note  that  Figure  3. IS  was  constructed  by  the  approximation  of  potential  vorticity.  q 

~    f  dT'dz.   whereas   Figure   3.17  was   constructed   from  the   full   form  of  potential 

vorticity   (eqn.    3.3).     The    neglected    term   in    Watts'    approximation   was    found    to 

contribute   between    10   and   20%    of  the   potential   vorticity  in   most   cases   and.   in 

particular,  the  second  term  dominated  in  experiment  1  on  the  near-shore,  cyclonic  side 
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Figure  3.17     Potential  vorticity  in  the  CCS  (exp.  1). 
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Figure  3.18     Same  as  Fisure  3.17  except  for  the  Gulf  Stream 
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of  the  coastal  jet.  The  scaled  form  of  potential  vorticity  was  used  in  the  Gulf  Stream 
analysis  due  to  limitations  in  data.  As  a  result  of  this,  only  qualitative  comparisons 
have  been  made  between  analyses. 

As  previously  discussed,  a  necessary  condition  for  instability  is  that  the  cross- 
stream  derivative  of  potential  vorticity  must  change  sign  somewhere  within  the  domain. 
To  examine  the  potential  vorticity  field  for  this  condition,  the  horizontal  gradient  was 
computed  (dq  dx)  and  multiplied  by  one  grid  length  (AX)  for  scaling  purposes.  The 
resultant  cross-stream  derivative  is  displayed  in  Figure  3.19  for  experiment  1.  As 
expected  from  Figure  3.17".  the  cross-stream  derivative  changes  sign  in  the  upper  levels 
across  the  extent  of  the  simulated  coastal  jet.  The  necessary  condition  for  instability  of 
the  mean  flow  was  therefore,  easily  satisfied  in  experiment  1. 

The  other  experiments  conducted  in  this  study  displayed  very  similar  patterns  of 
the  cross-stream  gradient  of  potential  vorticity  and.  therefore,  all  will  not  be  shown. 
All  of  the  experiments  satisfied  the  necessary  condition  for  instability.  There  was  no 
significant  difference  in  the  gradient  of  potential  vorticity  between  the  experiments  with 
uniform  wind  and  those  which  included  wind  stress  curl.  However,  the  gradient  of 
potential  vorticity  in  the  experiments  with  zero-slip  was  stronger  (both  positive  and 
negative)  than  those  with  the  free-slip  boundary  conditions  because  of  the  additional 
contribution  from  horizontal  shear  as  discussed  in  section  III.B.l.  (see  Figure  3.20). 

Why.  then,  if  the  necessary  condition  for  instability  has  been  satisfied  in  each 
experiment,  has  the  mean  flow  not  gone  unstable?  The  answer  to  that  question  is  not 
readily  apparent.  It  is  useful,  though,  to  explore  a  simple  baroclinic  instability  model 
to  determine  what  the  sufficient  conditions  are  for  baroclinic  instability  of  the  mean 
flow. 

The  simplest  model  which  can  incorporate  baroclinic  processes  is  a  two-layer 
model.  Following  the  method  of  Holton  (1979),  the  ocean  is  divided  into  two  discrete 
layers.  In  all  of  the  experiments  conducted  in  this  study,  the  southward-flowing  coastal 
jet  extended  to  approximately  1000  m  in  the  vertical.  For  this  reason,  the  upper  model 
layer  extends  from  the  surface  to  500  m  depth  and  the  lower  layer  extends  from  500  to 
1000  m.  See  Table  3  for  the  arrangement  of  variables  in  the  vertical  for  the  two-layer 
baroclinic  model. 

The  QG  vorticity  equation  is  applied  at  levels  1  and  3  whereas  the  QG 
thermodynamic  equation  is  applied  at  level  2.  The  streamfunction.  \\f~  is  obtained  by 
linearly  interpolating  between  -250  m  and  -750  m.  This  forms  a  closed  set  of  prediction 
equations  with  P  =  0,  as  in  the  PE  model. 
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Figure  3.19     Cross-stream  derivative  of  potential  vorticity 

experiment  1. 
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Figure  3.20     Same  as  Figure  3.19  except  experiment  3. 
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TABLE  3 

ARRANGEMENT  OF  VARIABLES  IN  THE  VERTICAL 
FOR  THE  TWO-LAYER  BAROCLIMC  MODEL 


0  m  wQ  =  0  level  0 

•250  m  \|/1  level  1 

-500  m  w,  level  2 

-750  m  \\J 3  level  3 

-1000  m  w4  =  0  level  4 


The  length  scale  used  in  this  model  is  the  Rossby  radius  of  deformation  (X*1) 
defined  by: 

X2  -  JL-.  ,  (3.5) 

d  (Az)2 

where  Az  is  the  layer  thickness  (500  m)  and  cr  is  the  stability  parameter  defined  by: 

<T  =  a  g  3T/dz  .  (3.6) 

The  phase  speed  (c)  is  obtained  from  solving  the  closed  set  of  QG  equations,  by 
assuming  wave-type  solutions,  yielding: 


k2  -  2X2     ,,  , 


1/2 


c       =        VM        ±        VT       [   _-_  ,J*i'  ,  (3.7) 


where  the  vertically  averaged  current  (Vu)  is  defined  as: 

V,  +  V 
M 


VVf        ■       11- 3    .  (3.S) 
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and  the  basic  state  "thermal  wind"  (Vj)  is  defined  as: 

V.  -  V, 
VT        =       -L— 1      .  (3.9) 

For  waves  with  meridional  wavenumbers  (k)  satisfying  k.  <  2X  ,  there  exists  an 
imaginary  component  of  the  phase  speed  (c).  This  criteria  defines  the  critical 
wavelength  (L  )  as  follows: 

L    =  21  2  (tt.X)   .  (3.10) 


or 


J/2 


L    =  [  (Az)  n  (2o)ilz  ] !  L   .  (3.11) 


Since  P  =  0,  the  criteria  for  instability  does  not  depend  on  the  magnitude  of  the  basic 
state  "thermal  wind"  (Vj).  All  wavelengths  longer  than  Lc  are  unstable  even  for  very 
small  vertical  shear. 

Any  wave  longer  than  this  critical  wavelength  will  amplify  at  the  growth  rate  (a) 
where 

a  =  kc    .  (3.12) 

The  e-folding  growth  time  is  simply  the  inverse  of  a.  Note  that  a  depends  on  Vy  and 
(T.  and  is  a  function  of  the  meridional  wavenumber,  k.  It  is  a  maximum  at  an 
intermediate  wavenumber  (0  <  k"  <  2K  ). 

Table  4  shows  the  inverse  Rossby  radius  of  deformation,  the  critical  wavelength 
(L  ),  the  basic  state  "thermal  wind",  and  the  most  unstable  wavelensth  (L  )  with  its 
associated  e-folding  time  for  experiments  1-4  and  for  an  imposed,  unstable  baroclinic 
jet  experiment  (BCJET)  from  a  previous  flat  bottom  study  (Batteen  et  al.,  1985).  In 
this  analysis,  the  computations  were  made  using  the  time  (days  81-90),  zonal  (within 
150  km  of  the  eastern  boundary),  and  meridionally  averaged  data  from  the  PE  model 
experiments. 

The  e-folding  times  for  all  of  the  experiments,  at  the  most  unstable  wavelengths, 
are  very  comparable.  The  computed  values  of  the  e-folding  times  are  only  an 
approximation  to  the  true  baroclinic  growth  rate  due  to  the  simple  two-layer  baroclinic 
instability  model  used,  and  therefore  the  difference  in  e-folding  time  between  all  of  the 
experiments  is  probably  negligible. 
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TABLE  4 

TWO-LAYER  BAROCLIMC  MODEL  RESULTS 

exp. 

\-[{km) 

L  c(km) 

Vj  (cm/s) 

L     (km) 

a"1  (days) 

#1 

16.96 

75.3 

12 

115.3 

2.79 

#2 

16.96 

75.3 

8.5 

115.3 

3.94 

17.14 

76.1 

11 

116.1 

3.08 

^4 

16.39 

72.8 

8 

112.8 

4.05 

BCJET 

10.9 

4S.4 

5.5 

78.4 

3.92 

The  short  e-foldmg  time,  however,  would  indicate  that  eddy  development  should 
be  apparent  after  about  a  week,  or  so  of  model  integration  .once  this  unstable  state  had 
been  achieved.  Indeed,  eddy  generation  was  apparent  in  the  experiment  with  the 
imposed  baroclinic  jet  after  about  25  days  of  model  integration.  Notice  that  there  is  a 
large  difference  in  the  critical  wavelength  and  associated  Rossby  radius  of  deformation 
between  experiments  1-4  and  the  imposed  baroclinic  jet  experiment.  Possibly  the 
reason  eddies  have  not  been  generated  in  the  wind-forced  experiments  lies  in  the 
difference  between  critical  wavelengths. 

Figure  3.21  is  a  plot  of  the  e-folding  times  as  a  function  of  meridional  wavelength 
(L  =  2tt  k)  for  the  imposed  baroclinic  jet  experiment  and  experiments  1-4.  The  two- 
layer  baroclinic  model  was  applied  over  a  ten-day  averaged  period  of  model  days  S I  -90. 
The  difference  in  critical  wavelengths  is  readily  apparent  in  the  figure.  Additionally, 
the  most  unstable  wavelength  is  ~  SO  km  in  the  imposed  baroclinic  jet  experiment 
versus  ~  115  km  in  the  windforced  experiments. 

To  determine  the  baroclmicity  of  the  experiments  throughout  the  spin-up  period, 
ten-day  averages  were  computed  and  the  two-level  baroclinic  model  was  applied  for 
three  different  periods  in  experiment  1.  Figure  3.22  shows  this  progression  and  it  is 
clear  that  the  baroclinic  How  in  this  experiment  is  becoming  more  unstable  throughout 
model  spin-up.  At  days  21-30  the  shorter  wavelengths  had  a  definite  preferred  growth 
rate  while  at  davs  81-90,  the  longer  wavelengths  are  not  considerably  less  unstable  than 
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Figure  3.21     Baroclinic  instability  growth  rates 
tor  all  experiments." 
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the  shorter  wavelengths.  Superimposed  on  Figure  3.22  is  the  familiar  result  from  the 
imposed  baroclinic  jet  experiment.  This  figure  seems  to  imply  that  if  the  experiments 
were  continued  eddies  may  be  generated.  The  problem  with  model  continuation  is  that 
the  desired  realism  would  be  lost.  Re-design  of  model  spin-up  would  be  required 
including  time  variant  wind  forcing  and  variable  parameterization  of  surface  heat  fiuxes 
to  simulate  conditions  indicative  of  the  fall  season.  Possibly,  the  best  solution  would 
be  to  simulate  an  annual  cycle  to  spin  the  model  up.  Additional  computer  resources 
would  be  required  for  an  experiment  of  that  scale. 

Once  a  stability  analysis  of  the  mean  flow  is  completed,  it  is  necessary  to  examine 
model  heat  and  momentum  diffusion  and  its  associated  damping  (e-folding)  time.  If 
the  damping  time  is  on  the  order  of  the  growth  rate  for  the  unstable  wavelengths, 
instability  of  the  mean  flow  will  be  suppressed. 

The  PE  model  used  in  this  study  incorporates  biharmonic  heat  and  momentum 
diffusion.  Following  Holland  and  Batteen  (1986),  the  QG  baroclinic  mode  vorticity 
and  thermal  equation  has  the  form: 

(V2-'k2)\\Jl    =    ... -A  V6y    4-    X2  B  V4  y    ,  (3.13) 

where  \j/  is  QG  "temperature"  (vj/,-V3  in  two-layer  model),  A  is  biharmonic  eddy 
viscosity,  B  is  biharmonic  eddy  diffusion,  X  is  inverse  Rossby  radius  (as  before),  and 
only  the  damping  terms  are  kept  in  (3.13). 

Assuming  wave  numbers  (k,  1),  (3.13)  becomes: 

A(k2  +  I2)3  +  X2B(k2  +  l2)2 
*.    -     "    I  k;  +  ,:  ;  x:  1    *  .  (3.14) 

or  \|/  =  \\fQ  e"Yl   .  (3.15) 

The  damping  (e-folding)  time  of  baroclinic  modes  is  therefore 

,  (k2  +  l2  4-  X2) 

Y"1      =  ,        ,  , 5 - ,  ,  .  (3.16) 

;  A(k2  +  l2)3  4-  X2B(k2  +    l2)2 

Table  5  shows  damping  rates  (y~  )  and  time  scales  y  for  various  wavelengths 
using  values  of  A  =  B  =  2  x  10  cm4  s"1  as  in  the  PE  model.  It  can  be  seen  that 
for  the  most  unstable  wavelengths  predicted  by  the  two-layer  baroclinic  model  for 
experiments    1-4   (~    100   km),    the   damping   time   scale   is    ~    95    days.     Since   the 
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baroclinic  instability  e-folding  time  was  previously  shown  to  be  ~  3-4  days,  diffusive 
damping  would  appear  to  be  negligible  at  that  wavelength  compared  with  the 
baroclinic  growth  rate.  For  longer  wavelengths,  which  also  have  small  e-folding  times 
due  to  baroclinic  instability,  the  diffusive  damping  is  entirely  negligible. 


TABLE  5 

DAMPING  RATES  FOR  VARIOUS  WAVELENGTHS 
IN  THE  PE  MODEL 

L  (km) 

Y  (days'1) 

Y"1  (days) 

10 

1.11    x    102 

.009 

20 

0.67   x    101 

.15 

40 

4.17   x    10'1 

2.4 

80 

2.63   x    io  "2 

3S 

100 

1.0S   x    uy2 

93 

200 

6.76  x    10"4 

14S0 
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IV.  COMPARISON  OF  MODEL  RESULTS  WITH  OBSERVATIONS 

A.       HISTORICAL  NOTES  ON  THE  CALIFORNIA  CURRENT  SYSTEM 

A  review  and  synthesis  of  existing  literature  on  the  CCS  through  1978  is 
presented  by  Hickey  (1979).  Throughout  the  history  of  observations  in  the  CCS,  a  set 
of  names  have  been  assigned  to  the  various  poleward  and  equatorward  flows  of  the 
CCS.  The  California  Current  is  the  equatorward  flow;  the  California  Undercurrent 
refers  to  poleward  flow  over  the  slope;  the  Davidson  Current  is  the  poleward  flow 
north  of  Point  Conception  during  fall  and  winter;  and  the  Southern  California 
Countercurrent  refers  to  the  poleward  flow  south  of  Point  Conception  in  the  California 
Bight  (Hickey,  1979).  Hickey  (1979)  points  out  that  as  our  sampling  resolution 
increases,  confusion  in  terminology  arises  and  that:  "The  confusion  in  nomenclature 
parallels  a  growing  uncertainty  about  the  temporal  and  spatial  variability  of  the  various 
currents  and  the  dynamical  relationships  among  them."  The  decade  following  197S  has 
yet  to  provide  clear-cut  definitions  of  the  structure  of  the  CCS  and  has  only  helped  to 
complicate  the  picture  by  demonstrating  that  the  CCS  is  more  variable,  both  spatially 
and  temporally,  than  previously  thought. 

Regardless  of  the  complex  variability,  the  flow  which  has  been  denoted  the 
California  Current  generally  has  both  an  offshore  and  a  nearshore  maximum  in  its 
alongshore  component.  The  large  scale  equatorward  flow  is  believed  to  be  fed  by  the 
northern  limb  of  the  North  Pacific  gyre  (the  West  Wind  Drift).  This  large  scale 
alongshore  component  of  southward  flow  has  been  observed  to  exceed  distances  of  500 
km  offshore  in  some  locations.  The  nearshore  maximum,  or  coastal  jet,  escaped 
attention  in  early  studies  of  the  CCS  primarily  because  the  majority  of  hydrographic 
data  had  station  spacing  of  about  50  km.  Recently,  however,  the  emphasis  has  shifted 
to  conducting  surveys  of  finer  resolution  in  an  attempt  to  resolve  the  synoptic- 
mesoscale  structure  of  the  coastal  jet  and  its  associated  meanders,  seaward  jets,  and 
eddies. 

In  this  study  of  wind  forcing  experiments  in  the  CCS,  the  large-scale  component 
of  equatorward  flow  associated  with  the  North  Pacific  Gyre  has  not  been  included. 
The  experiments  were  conducted  with  an  open  boundary  regional  model  with  local 
wind  forcing  and  for  clarity  of  interpretation,  large  scale  gyre  circulation  was  neglected. 
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In  future  developments  of  this  model,  it  may  become  necessary  to  include  large  scale 
gyre  circulation  in  the  quest  for  realistic  simulations  of  the  CCS.  For  the  above 
reasons,  only  observations  of  the  coastal  jet  and  its  associated  undercurrent  will  be 
used  for  comparison  with  model  results. 

B.       THE  COASTAL  JET 

In  Chapter  III.  the  model  spin-up  phase  was  described  in  detail.  That 
description  was  based  on  classical,  two-dimensional  coastal  upwelling  processes.  The 
resultant  coastal  jet  has  been  a  prominent  feature  in  descriptions  of  sheif  circulation  in 
the  Pacific  Northwest  during  the  upwelling  season  (Hickey.  1979).  Actual  descriptions 
of  the  coastal  jet,  however,  indicates  that  it  varies  widely  in  space  and  time.  Recently, 
hydrographic  and  X3T  surveys  have  demonstrated  that  instantaneous  cross-sections  of 
the  California  coastal  region  show  great  variability.  Kosro  (1986)  concluded  that 
synoptic  maps  of  near-surface  currents  off  northern  California  deviated  substantially 
from  the  classical  description  of  two-dimensional  wind  driven  upwelling.  Huyer  and 
Kosro  (1986)  conducted  six  surveys  near  Point  Arena.  California  during  the  upwelling 
season  in  19S1  and  1982.  The  simple  two-dimensional,  equatorward  baroclinic  coastal 
jet  was  observed  during  only  one  of  the  six  surveys  (20-22  April  1982).  The  simple 
two-dimensional  coastal  jet  is  almost  invariably  resolved  in  the  averaged  velocity  fields 
off  California  during  the  upwelling  season.  Huyer  and  Kosro  (1986)  found  that  the 
average  velocity  fields  (computed  over  all  6  surveys)  consistently  showed  baroclinic, 
equatorward  flow  in  the  upper  layers  and  weak  poleward  flow  along  the  shelf-break  at 
a  depth  of  ~  150  m.  Additionally,  it  was  noted  that  the  eddies  and  meanders  which 
dominated  the  velocity  field  during  most  surveys  were  absent  from  the  average  field. 
All  of  their  surveys  revealed  a  baroclinic  equatorward  jet  which  wound  around  eddies 
and  at  times  turned  directly  offshore.  It  is  hypothesized  that  since  the  classical  form  of 
the  coastal  jet  appears  in  the  average  velocity  fields,  but  only  rarely  in  a  synoptic 
survey,  the  upwelling  induced  coastal  jet  becomes  unstable  off  northern  California 
during  the  upwelling  season. 

The  classical  two-dimensional  coastal  jet  is  somewhat  more  prevalent  in  its  pure 
form  off  the  coast  of  Oregon.  In  two  successive  years  (1965  and  1966),  the  frontal 
zone  of  the  coastal  upwelling  region  off  Oregon  was  studied  by  Mooers  et  al.  (1976). 
They  found  that  the  computed  velocity  field  showed  equatorward  surface  flow  and  a 
poleward  undercurrent  below  the  inclined  frontal  layer.    Both  the  southward  surface 


flow  and  the  northward  undercurrent  exhibited  a  "jet-like  structure."  The  southward 
jet  was  most  intense  (  ~  20  cm  s)  within  20  km  from  the  shore  and  in  the  upper  40  m. 
The  poleward  undercurrent  was  most  intense  (  ~  10  cms)  between  10  and  40  km 
offshore  and  at  a  depth  of  SO  to  200  m. 

Huyer  and  Kosro  (1986)  state  that,  as  previously  discussed,  a  baroclinic  coastal 
jet  was  manifested  in  classical  form  during  their  20-22  April  19S2  survey  oil  northern 
California.  The  strength  and  position  of  the  jet's  core  (over  the  shelf-break)  were 
nearly  uniform  over  the  entire  alongshore  extent  of  the  survey  region  (~  100  km). 
This  anisotropic  characteristic  (more  rapid  variation  cross-shore  than  alongshore)  is  a 
recurrent  feature  of  the  coastal  jet  in  its  stable  form.  The  coastal  jet  was  further 
characterized  by  a  clearly  defined  core  located  about  25  km  from  the  coast  with 
southward  velocities  >  50  cm,  s  at  20  m  depth.  Increasing  in  depth  to  SO  m.  the  core 
was  located  farther  from  shore  ( ~  40  km)  with  maximum  southward  velocities  less 
than  30  m/s.  Available  data  suggested  that  there  existed  a  narrow  (10-20  km)  poleward 
undercurrent  at  ~  150  m  depth  with  maximum  velocities  of  about  15  cm/s. 

Strong  horizontal  and  vertical  shear  is  a  persistent  characteristic  of  the  baroclinic 
coastal  jet.  Huyer  (1983)  found  the  maximum  vertical  shear  to  be  relatively  constant 
(~  4  x  10°  s"1)  during  the  upwelling  season  in  spite  of  significant  wind  stress 
fluctuations.  Huyer  and  Kosro  (19S6)  used  current  meter  data  to  determine  that 
average  maximum  vertical  shear  was  4-5  x  10  s"1  in  the  upper  35  m  over  the  outer 
shelf  during  the  19S1  and  1982  surveys.  Kosro  (1986)  computed  vertical  shear  to  be 
greater  than  2  x  10  s"1  in  the  coastal  jet  off  northern  California  during  CODE  with 
observed  cyclonic  horizontal  shear  exceeding  10  s"1  which  is  the  local  value  of  the 
Coriolis  parameter  (f).  Mooers  et  al.  (1976)  computed  mean  geostrophic  vertical  shear 
from  the  thermal  wind  equation  to  be  ~~  4  x  10"3  s"1  ofT  Oregon  which  compared  well 
with  their  measured  current  shear  of  4.6  x  10"3  s"1  in  the  coastal  jet  off  Oregon. 
Kosro  and  Huyer  (1986)  collected  data  across  a  seaward  jet  off  northern  California  and 
suggested  that  it  was  continuous  with  the  coastal  jet.  Horizontal  shear  on  the  cyclonic 
side  of  the  jet  was  estimated  to  be  on  the  order  of  10  s"1  (f)  and  on  the  order  off/3  or 
less  for  the  anticyclonic  side  of  the  jet. 

C.       COMPARISONS 

The  most  important  consideration  in  a  comparison  oi'  model  results  with 
observations  is  which  model  assumptions  might  cause  the  model  results  to  deviate  most 
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significantly  from  nature.  In  this  study  there  are  two  rather  artificial  assumptions,  the 
flat-bottom  and  the  straight  coastline,  which  could  significantly  alter  the  results  from 
that  observed  in  nature.  As  discussed  in  Chapter  II.  however,  the  flat-bottom  and 
straight  coastline  are  crucial  in  interpretation  of  early  modeling  results.  It  is  with  this 
in  mind  that  the  following  comparisons  are  presented.  Table  6,  below,  is  a  synopsis  of 
some  of  the  more  interesting  features  of  the  coastal  jet  generated  in  experiments  1-4 
along  with  generalized  features  of  a  classical,  upwelling-induced,  two-dimensional 
coastal  jet  occasionally  found  off  the  coast  of  northern  California  and  Oreson. 


TABLE  6 

STATISTICAL  COMPARISON  OF  MODEL  RESULTS 
WITH  OBSERVATIONS  OF  THE  COASTAL  JET 

exp.l 

exp.2 

exp.3 

exp.4 

obs. 

A. 

100  cm/s 

75  cm/s 

100  cm/s 

70  cm/s 

30-50  cm/s 

B. 

25  km 

25  km 

25  km 

25  km 

~  25  km 

C. 

20  km 

20  km 

28  km 

28  km 

"-  35  km 

D. 

100  km 

190  km 

100  km 

190  km 

25-100  km 

E. 

1000  m 

1000  m 

1000  m 

1000  m 

<  100  m 

F. 

5-10  cms 

5-10  cm 

5    5-10  cm/s 

5-10  cm  s 

5-15  cm  s 

G. 

60  km 

25  km 

50  km 

20  km 

~  15  km 

H. 

2.8  s'1 

2.6  s"1 

2.9  s"1 

2.6  s_1 

2-5  s_1 

I. 

,15f 

.If 

.6f 

,4f 

~f?2 

J. 

-,2f 

-.15f 

-.2f 

-.15f 

~-f?3 

K. 

.07  °/km 

.05  7km 

.06  °/km 

.05  7km 

-.05  7km 

L. 

8.5  °C 

9.5  °C 

9.0  °C 

10.0  °C 

12-13  °C 
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The  letter  designations  in  Table  6  have  the  following  meanings: 

A.  Maximum  jet  velocity. 

B.  Offshore  location  of  surface  jet  axis. 

C.  Offshore  location  of  depth  averaged  jet  axis. 

D.  Offshore  extent  of  jet. 

E.  Depth  of  jet. 

F.  Maximum  undercurrent  velocity. 

G.  Width  of  undercurrent. 

H.    Maximum  vertical  shear  (  x    10"-). 
I.    Maximum  horizontal  cyclonic  shear. 
J.    Maximum  horizontal  anticyclonic  shear. 
K.    Xearshore  surface  temperature  gradient. 
L.    Mean  coastal  surface  temperature. 

The  maximum  jet  velocities  that  developed  in  experiments  1-4  exceed  the 
observed  values  by  about  a  factor  of  two.  The  exceptionally  strong  surface  currents 
developed  in  experiments  1  and  3  may  be  explained  in  part  by  the  unrealistically  large 
value  of  wind  stress  at  the  coast  in  those  experiments.  The  surface  stress  values  were 
chosen  by  taking  an  offshore  mean  climatological  value  and  applying  it  uniformly 
across  the  model  domain.  Experiments  2  and  4,  which  included  the  curl  of  wind  stress, 
had  more  realistic  (reduced)  values  of  wind  stress  near  the  coast  and  indeed  produced 
more  realistic  surface  currents.  Additionally,  winds  were  input  as  steady  which  may 
produce  stronger  currents  than  a  time  varying  forcing  function.  Time  variable  winds 
may  by  a  key  to  generating  instabilities  in  the  coastal  jet  and  thus  limiting  its  growth. 
Note  that  experiment  4  produced  the  maximum  jet  velocities  closest  to  those  observed. 

The  offshore  location  of  the  jet  axis,  at  the  surface  and  with  depth,  compared 
favorably  with  observed  locations.  The  vertical  structure  of  the  coastal  jet  tends  to 
follow  the  coastal  slope.  The  vertical  eastern  boundary  in  the  model  produced  an 
almost  vertically  stacked  jet,  whereas  the  slope  off  northern  California  produces  a 
coastal  jet  which  tilts  offshore  with  depth.  Associated  with  the  axis  location  is  the 
offshore  extent  of  the  jet.  Observations  indicate  that  the  coastal  jet  is  a  narrow  feature 
which  extends  horizontally  offshore  from  25-100  km.  Experimental  results  show  a 
coastal  jet  about  100  km  wide.  The  inclusion  of  wind  stress  curl  (experiments  2  and  4) 
tended  to  increase  the  width  of  the  modeled  equatorward  flow  but  the  main  body  of 
the  jet  remained  about  100  km  in  width. 
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The  depth  of  jet  penetration  in  the  water  column  is  by  far  the  largest  deviation  of 
model  results  from  observations.  Further  experimentation  and  analysis  is  needed  to 
identify  the  cause  of  this  major  discrepancy.  One  possibility  is  the  absence  of  a 
shelf  slope  topography  in  the  model.  Model  basin  depth  is  constant  at  4500  m  which 
means  that  the  coastal  jet  penetrated  to  about  one-fourth  of  the  basin  depth. 
Observations  show  that  the  observed  coastal  jet  penetrates  to  about  one-third  of  the 
continental  shelf  depth.  Therefore  a  relative  comparison  of  model  results  with 
observations  is  good.  Preliminary  results  of  a  wind  forcing  experiment  with  idealized 
cross-shore  topography  indicates  that  the  model  coastal  jet  penetration  is  in  fact 
dependent  on  shelf  depth.  An  additional  consideration  in  the  vertical  structure  of  the 
jet  is  the  vertical  temperature  stratification.  McCreary  et  al.  ( 1 9 S 6 )  suggests  that  a 
shallow  (~  100  m)  seasonal  thermocline  may  be  very  important  in  the  formation  of 
the  coastal  jet,  whereas  the  model's  basic  stratification  was  that  corresponding  to  the 
major  thermocline  having  a  vertical  scale  of  450  m. 

The  poleward  undercurrent  maximum  velocities  in  experiments  1-4  were 
somewhat  slow  as  compared  with  observations  and  also  seem  independent  of  the 
magnitude  of  wind  stress.  Topography,  again,  may  be  an  important  mechanism  in  the 
formation  of  undercurrents.  Open  boundary  conditions  in  the  model  may  modify  the 
formation  of  an  undercurrent  as  compared  with  the  effects  of  gyre  scale  circulation  in 
the  North  Pacific.  Wind  relaxation  events  off  the  coast  of  California  have  also  been 
considered  to  be  closely  associated  with  undercurrents  and  countercurrents  within  the 
CCS.  Mooers  et  al.  (1976)  stated  that  "The  relationship  between  the  northward 
Davidson  Inshore  Current,  the  California  Countercurrent,  and  the  undercurrent 
associated  with  coastal  upwelling  is  not  understood."  It  is  probably  safe  to  say  that 
that  statement  still  holds  true  today.  Recent  emphasis  in  the  oceanographic 
community  has  been  placed  on  efforts  to  understand  those  relationships.  The  width  of 
the  undercurrent  is  generally  considered  very  narrow  (~  15  km),  as  described  by  Huyer 
and  Kosro  (19S6).  Experiment  4  produced  the  narrowest  undercurrent  (20  km)  which 
corresponds  well  with  observations. 

Vertical  and  horizontal  shear  measurements  have  demonstrated  that  the  coastal 

jet  off  California  and  Oregon  is  a  highly  sheared  system.    Values  for  the  vertical  shear 

range  from  about  2   x    icr3  s"1  to  5    x    10"3  s"1  (Huyer,  19S3;  Kosro,  19S6:  Mooers  et 

al.,  19S6).    Experimental  shear  values  fall  within  that  range  but  on  the  low  side  (~   3 

x    10     s"1).    The  topographic  influence  may  again  be  important  to  the  establishment 
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of  large  vertical  shear.  The  coastal  shelf  and  slope,  in  compressing  the  vertical 
structure  of  the  jet.  may  also  increase  the  associated  vertical  shear.  The  measurements 
of  horizontal  shear  in  the  CCS  suggest  that  a  typical  value  for  the  cyclonic  shear  may 
be  on  the  order  of  f/2  to  f  and  that  the  anticyclonic  shear  is  on  the  order  of  -f  3 
(Kosro,  1986).  Experimental  results  demonstrate  that  a  free-slip  boundary  condition 
on  the  eastern  boundary  (experiments  1  and  2)  prevents  the  development  of  large 
horizontal  cyclonic  shear.  Even  with  the  zero-slip  boundary  condition  (experiments  3 
and  4)  the  largest  cyclonic  shear  developed  by  the  PE  model  was  .6f.  Anticyclonic 
shear  was  similarly  smaller  in  the  model  experiments  than  the  observed  values  of 
~-f;3. 

Of  course,  the  coastal  temperature  structure  is  closely  tied  to  the  coastal  jet 
structure  owing  to  its  baroclinic  nature.  The  observed  offshore  surface  temperature 
gradient  during  the  upwelling  season  off  northern  California  is  on  the  order  of  .05  °  km 
with  a  mean  July  coastal  temperature  of  10-13  °C  (Huyer  and  Kosro,  19S6V  Model 
results  show  a  range  of  gradients  of  .05-. 07  °  km  with  coastal  temperatures  in  the  range 
of  S. 5-10  °C.  Those  values  compare  well,  and  the  cooler  model  coastal  temperatures 
may  be  explained  by  the  unrealistically  high  model  jet  velocities  or  inexact  surface  heat 
flux  parameterization. 

It  is  difficult  to  make  a  quantitative  comparison  of  experimental  results  with 
observations  of  the  coastal  jet  in  the  CCS.  Observed  values  of  jet  characteristics  vary 
widely  from  study  to  study  and  the  stable  form  of  the  classical  two-dimensional  coastal 
jet  has  been  infrequently  surveyed.  From  the  above  discussion,  it  is  clear  that  the 
structure  of  the  simulated  coastal  jet  is  sensitive  to  the  forcing  function,  boundary 
conditions,  heating  parameterization,  topography,  and  a  host  of  other  conditions  too 
numerous  and  complex  for  the  scope  of  this  study.  Experiment  4,  which  included 
idealized,  climatological  wind  stress  curl  with  a  zero-slip  boundary  condition  on  the 
eastern  boundary,  was  the  best  simulation  of  the  coastal  jet  of  the  four  experiments 
conducted  in  this  studv. 


77 


V.  SUMMARY  AND  CONCLUSIONS 

A.       SUMMARY 

This  study  has  examined  the  response  of  an  idealized,  flat-bottom,  eastern 
boundary  oceanic  regime  to  steady,  equatorward,  local  wind-forcing  on  a  seasonal  time 
scale.  Both  uniform  and  offshore-varying  wind,  with  zero-slip  and  free-slip  boundary 
conditions  on  the  eastern  boundary,  were  employed  as  forcing  mechanisms  in  a  ten- 
level,  high  resolution,  primitive  equation  ocean  model. 

A  classical  two-dimensional,  upwelling-induced,  equatorward  flowing  coastal  jet. 
with  an  associated  poleward  undercurrent,  developed  in  response  to  both  uniform  and 
offshore-varying  wind  forcing.  The  variations  in  formation  and  structure  of  the  coastal 
jet  and  undercurrent  were  explored  and  analyzed  in  four  numerical  experiments: 
uniform  wind  with  free-slip;  offshore-varying  wind  with  free-slip;  uniform  wind  with 
zero-slip;  and  offshore-varying  wind  with  zero-slip. 

1.  Stability  Analysis 

Instability  of  the  mean  flow  did  not  occur  in  any  of  the  aforementioned 
experiments  and.  therefore,  eddies  did  not  develop.  A  stability  analysis  was  conducted 
on  all  of  the  modeled  coastal  jets  in  this  study  along  with  an  imposed  coastal  jet  from  a 
previous,  unforced,  flat-bottom  experiment  (Batteen  er  ai,  1985)  which  formed  eddies 
after  approximately  thirty  days  of  model  run  time.  Potential  vorticity  was  computed 
throughout  the  model  domain  of  each  experiment  and  the  cross-stream  derivative  of 
potential  vorticity  was  found  to  change  sign  in  the  region  of  each  modeled  coastal  jet. 
This  fact  established  that  the  modeled  coastal  jets  all  satisfied  the  necessary  condition 
for  mixed  (barotropic  and  baroclinic)  instability. 

To  further  explore  conditions  for  instability,  a  two-layer,  quasigeostrophic. 
baroclinic  instability  model  was  applied  to  each  experiment.  It  was  established  that  all 
of  the  experiments  satisfied  the  sufficient  conditions  for  baroclinic  instability  as  defined 
by  the  simple  two-layer  stability  model.  The  four  experiments  conducted  in  this  study 
did  not  significantly  vary  in  critical  wavelengths,  most-unstable  wavelengths,  or  in 
baroclinic  instability  e-folding  times.  However,  the  imposed  baroclinic  jet  experiment 
(Batteen  et  ai,  1985)  had  a  shorter  critical  wavelength  and  a  shorter  most-unstable 
wavelength  than  the   four  wind-forced  experiments.     The  e-folding   rates   for  longer 
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wavelengths  were  somewhat  larger  in  the  imposed  baroclinic  jet  experiment  as 
compared  with  this  study's  experiments;  which  would  indicate  that  the  wind-forced  jets 
are  more  unstable  than  the  imposed  jet  at  longer  wavelengths.  The  reason  why 
instability  of  the  mean  How  did  not  occur  in  the  model  generated  coastal  jets  remains 
an  open  scientific  question. 

2.  Comparison  with  Observations 

A  comparison  of  model  results  with  observations  of  the  CCS  was  conducted. 
A  numerical  model  can  properly  simulate  nature  only  to  the  extent  which  observations 
have  provided  the  modeler  with  understanding  of  the  feature  he  desires  to  simulate. 
This  is  a  problematic  area  in  modeling  the  CCS.  Observations  of  the  CCS  demonstrate 
great  variability  in  time  and  space  and  the  classical  two-dimensional  upwelling-induced 
coastal  jet  has  been  infrequently  observed  and  studied  off  California.  With  this  in 
mind,  a  synthesis  of  literature  addressing  the  coastal  jet  was  conducted  and  a  range  of 
values  for  some  of  the  more  important  characteristics  was  compiled  and  compared  with 
model  results. 

The  model  generated  coastal  jet  compared  favorably  with  observations  in 
offshore  axis  location,  undercurrent  strength  and  width,  vertical  and  horizontal  shears 
and  surface  temperature  gradients.  Some  deficiencies  were  noted.  These  included 
discrepancies  in  the  strength,  depth,  and  offshore  extent  of  the  modeled  coastal  jet 
compared  with  observations.  Experiment  4.  which  included  wind  stress  curl  with  a 
zero-slip  boundary  condition,  is  believed  to  be  the  most  realistic  model  jet  produced  in 
this  study. 

3.  Effects  of  Offshore  Varying  Wind  Forcing 

Positive  wind  stress  curl,  within  200  km  of  the  eastern  boundary,  was  applied 
by  the  offshore-varying  wind  forcing  in  two  of  the  experiments  (experiments  2  and  4). 
The  positive  wind  stress  curl  induced  cyclonic  relative  vorticity  to  the  system  and  added 
an  additional  component  of  upwelling,  due  to  the  effect  of  Ekman  pumping,  to  the 
coastal  upwelling  component  (caused  by  equatorward  wind  forcing  along  the  eastern 
model  boundary).  The  Ekman  pumping  induced  upward  vertical  motion,  from  surface 
divergence,  within  200  km  of  the  coast  and  created  an  extended  zone  of  baroclinicity  in 
the  upper  layers  of  the  ocean.  The  extended  zone  of  baroclinicity,  therefore,  increased 
the  width  of  equatorward  flow  over  that  generated  by  the  experiments  with  uniform 
wind  forcing  (experiments  1  and  3).  Positive  wind  stress  curl  near  the  coast  also  moved 
the  depth-averaged  jet  axis,  along  with  the  undercurrent  axis,  farther  offshore  than 
forcing  without  curl. 
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4.  Effects  of  the  Zero-Slip  Boundary  Condition 

As  previously  discussed,  both  free-slip  and  zero-slip  boundary  conditions  were 
employed  on  the  eastern  model  boundary.  The  free-slip  boundary  condition  implies 
that  the  offshore  gradient  of  tangential  velocity  vanishes  at  the  boundary  (experiments 
1  and  2),  and  the  zero-slip  condition  sets  the  tangential  velocity  itself  to  zero  at  the 
coast  (experiments  3  and  4).  The  zero-slip  boundary  condition,  therefore,  imparted 
cyclonic  relative  vorticity  to  the  inshore  side  of  the  model  coastal  jet.  From  the 
comparison  of  model  results  with  observations  of  the  CCS.  it  is  believed  that  eastern 
boundary  coastal  circulation  is  best  simulated,  in  a  high  resolution  model,  with  the  use 
of  a  zero-slip  boundary  condition  on  the  eastern  boundary. 

B.       RECOMMENDATIONS 

This  study  was  a  first  step  in  a  numerical  modeling  effort  at  the  Naval 
Postgraduate  School  in  implementing  wind  forcing  into  a  PE  model  designed  to 
simulate  the  CCS  off  northern  California.  It  is  important  for  future  studies  to 
determine  the  reasons  why  instability  of  the  mean  flow  produced  by  wind-forcing,  in 
this  study,  does  not  occur  when  stability  analyses  indicates  that  the  flow  is  unstable. 
The  development  and  evolution  of  eddies  observed  in  the  CCS  is  of  extreme  interest, 
and  understanding  instability  of  the  mean  flow  may  be  the  key  to  understanding  eddy 
processes.  Designing  future  experiments  to  understand  the  vertical  structure  of  the  jet 
and  why  the  model  jet  is  much  stronger  and  wider  than  the  observed  coastal  jet  should 
be  given  high  priority  as  well. 

This  study  had  the  constraint  of  a  flat-bottom  basin  with  a  regular  coastline. 
Implementation  of  topography  first  requires  the  addition  of  depth-averaged  flow  within 
the  PE  model  for  scientific  validity.  Once  the  depth-averaged  flow  is  incorporated,  the 
door  will  be  opened  for  a  series  of  studies  of  wind-forced  circulation  over  idealized 
topography.  The  influence  of  topography  is  thought  to  have  a  great  impact  on  coastal 
circulation.  Additionally,  an  irregular  coastline  which  would  simulate  capes  and  bays 
along  the  coast  of  California  would  be  a  fascinating  study. 

Finally,  one  should  consider  the  inclusion  of  gyre-scale  circulation  as  the  baseline 
for  wind-forcing  experiments  in  the  CCS.  The  gyre-scale  circulation  could  be  used  to 
initialize  the  model  prior  to  model  spin-up  by  local  wind-forcing.  Time-dependent 
winds  could  eventually  be  investigated  as  the  prelude  to  model  forcing  by  actual, 
observed  wind  records. 
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APPENDIX 
ANALYTIC  APPROXIMATION  OF  UPWELLING  TIME  SCALE 

The  time  scale  over  which  upvvelling  occurs  may  be  analytically  approximated 
through  use  of  the  principle  of  conservation  of  mass,  scale  analysis,  and  Ekman 
transport.  The  upwelling  rate  can  then  be  compared  with  the  PE  model  results  for 
verification  of  model  physics.  The  approximation  involves  no  diabatic  heating  and 
neglects  the  component  of  velocity  due  to  pressure  gradient  force.  A  control 
experiment  was  conducted  which  included  no  diabatic  processes  (surface  thermal 
damping)  and  which  had  uniform,  steady  wind  forcing.  That  experiment  will  be  used 
for  comparison  with  the  following  approximation.  All  symbols,  in  the  following 
derivation,  are  as  defined  in  Chapter  II  unless  otherwise  stated. 

To  arrive  at  a  scale  for  the  vertical  component  of  velocity,  the  equation  of 
continuity,  with  the  assumption  of  incompressibility,  is  written  as 

du  dv  d\v 

+     +     =     0      .  (A.l) 

dx  dy  dz 

Assuming  no  longshore  variation  in  v.  then  equation  (A.l)  may  be  reduced  to 

<3u  d\v 

—  +     —      =       0      .  (A.l) 
ox              oz 

Scale  analysis  of  equation  (A. 2)  yields  the  following  scale  for  w  (vertical  velocity): 

U  H 

W       ~       ,  (A. 3) 

where  UH  is  mass  transport  and  L    is  a  length  scale. 

Since  it  is  desired  to  isolate  the  effect  of  upwelling  caused  only  by  motion  of  the 
fluid  away  from  the  solid  boundary  driven  by  surface  stress,  Ekman  mass  transport 
(Mjr)  is  required  where 


MEx  ~  J-H°  Po  UE  dz 


x 


f 


which  may  be  rewritten  as 
i 


—      =  P.  u  H 
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or 

UH       =       -i-  .  (A.4) 

Using  equation  (A. 3).   W  is   substimied  into   equation  (A.4)  to   arrive   at   the 
following  relationship  between  vertical  velocity  and  surface  stress: 


W       ~  -f-  .  (AS) 

Po  {  L.x 

The  surface  stress  can  be  calculated,  through  use  of  the  bulk  aerodynamic  formula  with 
a  constant  drag  coefficient,  by 

*,      =      Pa  CD    V  2    • 

Substituting  this  expression  for  t  into  equation  (A. 5)  results  in  the  final  form  for 
vertical  velocity  (W)  of 

W       ~  "—77-  ■  (A.6) 

Po  l  L, 

Using  a  value  of  V  =  900  cm  s"1  and  L  =  14.2  km  (the  second  internal  Rossby  radius 
of  deformation  computed  by  the  method  of  Feliks  (1985)).  equation  (A.6)  yields  a  scale 
value  for  the  vertical  velocity,  generated  by  coastal  upwelling,  of 

W    ~    0.00947  cm/sec    =    S.18  m'day. 

Model  output  demonstrates  that,  in  the  upper  levels  of  the  ocean,  vertical 
velocity  decreases  by  80-90%  within  two  gridpoints  of  the  coast.  This  length  scale,  of 
16  km,  compares  well  with  the  second  internal  Rossby  radius  of  deformation.  The 
model  length  scale  of  vertical  velocity  does  not,  however,  compare  well  with  the  first 
Rossby  radius  of  deformation  which  was  computed  to  be  29.3  km.  The  horizontal 
length  scale  corresponding  to  the  second  Rossby  radius  of  deformation  is  appropriate 
because  the  model  vertical  velocity  has  a  maximum  close  to  the  sea-surface. 

A  comparison  of  the  rate  of  model  upwelling  with  that  obtained  from  an  Ekman 
balance  was  also  made  starting,  arbitrarily,  at  model  day  30  since  Ekman  balance  is  not 
well  established  during  the  initial  days  of  model  spin-up.  At  this  time,  the  depth  of  the 
9.03°  C  isotherm  (chosen  for  illustrative  purposes)  at  the  coast  was  169.4  m  (level  4). 


S2 


Using  the  vertical  velocity  approximation  computed  above,  the  9.03°  C  isotherm 
should  be  displaced  to  level  1  (12.9  m)  in  20.65  days.  Model  output  demonstrates  that, 
at  day  50.65,  coastal  temperature  at  level  1  is  9.19°  C  which  compares  quite  favorably 
with  the  analytically  approximated  temperature  of  9.03°  C. 
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